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ABSTRACT 



We show that the dependence of Lya absorption on environment leads to signifi- 
cant non-gravitational features in the redshift space power-spectrum of Lya selected 
galaxies. We derive a physically motivated fitting formula that can be included in 
clustering analyses, and use this to discuss the predicted features in the Lya galaxy 
power-spectrum based on detailed models in which Lya absorption is influenced by 
gas infall and/or by strong galactic outflows. We show that power-spectrum measure- 
ments could be used to study the astrophysics of the galaxy-IGM connection, and to 
measure the properties of outflows from star-forming galaxies. Applying the modified 
redshift space power-spectrum to a Lya survey with parameters corresponding to the 
planned Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), we find that 
the dependence of observed Lya flux on velocity gradient and ionising background 
may compromise the ability of Lya selected galaxy redshift surveys to constrain cos- 
mology using information from the full power-spectrum. This is because the effects 
of fluctuating ionizing background and velocity gradients effect the shape of the ob- 
served power-spectrum in ways that are similar to the shape of the primordial power- 
spectrum and redshift space distortions respectively. We use the Alcock-Paczynski test 
to show that without prior knowledge of the details of Lya absorption in the IGM, the 
precision of line-of-sight and transverse distance measurements for HETDEX will be 
1.3 — 1.7%, decreased by a factor of ~ 1.5 — 2 relative to the best case precision of 
0.8% available in a traditional galaxy redshift survey. We specify the precision with 
which modelling of Lya radiative transfer must be understood in order for HETDEX 
to achieve distance measurements that are better than 1%. 

Key words: cosmology: diffuse radiation, large scale structure, theory - galaxies: 
high redshift, inter-galactic medium 



1 INTRODUCTION 



The Lya emission line of galaxies provides a primary observ- 
able for discovering high redshift galaxies (e.g. Kashikawa 
et al l [20061 |Ouchi et aLl[20Tol |Iye et al.|[2b06, .Lehnertl 
et al.|20f0 |, for studying their starformation and inter-stellar 



medium (e.g. Verhamme et al. 2006 Dessauges-Zavadsky 



et al.|2010 Steidel et al.|20lT l, and for studying the ionisa 



tion state of the inter-galactic medium or IGM (e.g. 


Haiman 


& Spaans"1999||Malhotra & Rhoads|2004| |Kashikawa et al. 


2006; Dijkstra et al. 2007b 


\. In addition to measuring the lu- 



et al.||2006| [Kashikawa et al.,,2006, Ouchi ct al.^ 2008 



Cas- 



sata et al.[[2011[ [Ouchi et al.[[2010[ [Blanc et al.[[2010[ ), sam- 
ples have recently become large enough to enable studies of 
Lya galaxy clustering (e.g. Gawiser et al.|2007 Kovac et al 



[2007|[Orsi et al.|2008[[Guaita et al.|2010[ [Ouchi et al.|2010[ ). 
These clustering studies yield complementary information 
to the luminosity function, since clustering of galaxies can 
provide a direct estimate of the halo mass, independent of 
the life-time of star-burst activity. Comparison with the lu- 
minosity function therefore provides an avenue to estimate 
the overall efficiency and duration of star-formation activity 
(e.g. Nagamine et al.|2010 l. 



Unlike populations of Ly-break galaxies that are se- 
lected via broad band photometry, the observed brightness 
of a Lyof emitter is sensitive to its local extra-galactic en- 
vironment. In particular, Lya radiation that escapes from 
galaxies can be scattered out of the line-of-sight by neu- 
tral hydrogen atoms in inter-galactic medium (IGM) sur- 
rounding the galaxy. For example, prior to the completion 
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of reionisation the strength of the damping wing in the Lya 
absorption hne means that Lya galaxies should be more eas- 
ily detected inside the HIT regions that are thought to have 
been generated by clustered early star-forming galaxies (e.g. 
SantoS|,2004, jWyithe fc Loeb||2005| [Mesinger fc Furlanettol 
2008a I, leading to enhanced clustering that could provide a 
signature of patchy reionisation (IMcQuin n et al.|2007 Iliev 



et al.|2008| [Mesinger fc Furlanett o 2008b 

Following the conclusion of reionisation, the fraction of 
radiation that is scattered out of the line-of-sight is depen- 
dent on resonant absorption within the highly ionised IGM. 
In this regime, the fraction of Lya flux that is transmit- 
ted to the observer is dependent on quantities like the infall 
velocity, the local overdensity of mass, and the ionising back- 
ground. Any environmental dependence of observed flux at 
fixed intrinsic luminosity therefore leads to an environmen- 
tal dependence of the host halo mass of observed galaxies at 
fixed observed flux. This in turn leads to a dependence of 
the observed galaxy density on environment that differs to 
that expected from galaxy bias. Since clustering studies are 
performed in flux limited surveys, the net result is a modi- 
fication of the observed clustering of Lya selected galaxies 
( Zheng et al.|20lT l. This modification of the observed clus- 
tering is non-gravitational. 

Galaxies are thought to provide a (biased) tracer of the 
density field, and hence their clustering can be used to infer 
the statistical properties of the mass-density field on the 



large scales that probe cosmology (e.g. 


Cole et al. 


2005 


Eisenstein et al.||2005| 


Percival et al.||2007 Okumura et al.| 


2008 


Gaztanaga et al. 


2008 


Reid et al.|2010 


1. The gravita- 



tional contributions to this clustering can be studied analyti- 
cally, providing qualitative interpretation of the observations 



(e.g. Sheth et al. 20011. However the precision of modern 
galaxy redshift surveys has meant that N-body analyses are 
required to understand the observations in detail (e.g. 



Tin- 



ker et al.|2006, Eisenstein ot al.|2007||Seo et al.|2008[ ). While 
gravitational effects are thought to be the only mechanism 
influencing clustering at an observable level in traditional 
galaxy redshift surveys, studies of clustering among Lya se- 
lected galaxies will need to also include non-gravitational 
contributions. For example, the effect of peculiar velocity 
gradients on the observed clustering of Lya selected galax- 
ies will not be the same as in a usual galaxy redshift survey 



(Zheng et al. 20111. This is because both the observed lu- 
minosity and the observed redshift space density of the Lya 
galaxies depend on velocity gradient. 



Zheng et al. (20111 employed a numerical simulation 
and calculated the full radiative transfer of Lya photons 
to study the clustering of Lya selected galaxies aX z = 5.7. 
They argue that the effect of velocity gradients leads to line- 
of-sight clustering that is suppressed, in contrast to the en- 
hancement seen in traditional galaxy redshift surveys (e.g. 
Kaiser 1987 Peacock et al. 20011. In addition to identify- 



ing this new clustering effect, Zheng et al. (2011 1 argue that 



clustering will be enhanced transverse to the line-of-sight, 
leading to a clustering amplitude that is increased relative 
to expectations in the absence of Lya transmission effects. 

An important difference between the study of 'Zheng' 
et al. (20111 and previous related work (e.g. Santos 2004, 



Dijkstra et"ar]|2007a| |Orsi et al.|p008| |Dayal' et al.|j'2009 i 



is the implementation of full radiative transfer within their 
simulation ( Zheng et al.||2010 1 . Many previous studies have 



utilised a model in which the fraction of Lya photons that 
reach the observer is equal to exp(— r), where r is the in- 
tegrated Lya optical depth along the line of sight. However 



Zheng et al. (20111 argue that photons are scattered back 
into the line-of-sight along directions of low density, and so 
must be included in addition to the directly transmitted pho- 
tons accounted for in the exp(— r) model. They argue that 
while the exp(— r) model provides a qualitative explanation 
of observed Lya galaxy properties, quantitative differences 
are found which require full radiative transfer to interpret. 

The issue is partially related to the size, and hence the 
surface brightness, of the region from which Lya photons 
are scattered to the observer. This size is very sensitive 
to the properties of the gas within the virial radius of a 
galaxy, and also depends on the intrinsic broadness of the 
line ( Laursen et al.pOll l. In a recent paper Laursen et al. 
(20111 have modeled sight-lines to galaxies at lower redshift 
and at much higher resolution than are available in the sim- 



ulations of Zheng et al. (20101. In contrast to Zheng et al 



(20111, Laursen et al. (20111 flnd that the exp(— r) model 



provides an excellent description of the fraction of transmit- 
ted flux that they compute from their high resolution radia- 



tive transfer simulations. However Laursen et al. (20111 are 
not able to study the effect of transmission on Lya galaxy 
clustering. 

We predominantly consider clustering of Lya galaxies at 
z < 3. Following Laursen et al. (2011 1, we therefore employ 



an exp(— r) model for Lya transmission, and use this to ex- 
plore the possible effect of environmental dependence of Lya 
transmission on the clustering of Lya selected galaxies. This 
enables us to draw on the extensive prior work evaluating 
the possible effects of different astrophysical effects including 
infall and star-formation rates (Dijkstra et al. 2007a I, and 



galactic winds and outflows (Ahn et al. 2003 



Verhamme 



et al.|2008| [Dijkstra fc Wyithe|2010| ). The detailed analytic 



models we use in this work are in good agreement with the 



simulations of Laursen et al. (20111. Since we do not know 



the phenomena that are most important in setting the ob- 
served flux of Lya emitting galaxies, the goal of this work 
is not to produce a detailed model for the Lya luminosity 
function itself, or to predict the clustering amplitude (e.g. 
Orsi et al. 2008 Shimizu et al. 20111. Rather, we aim to 



understand the effect that fluctuations in transmission will 
have on the shape of the power-spectrum as a function of 
scale and direction, as well as on the clustering amplitude. 

Most studies of clustering among Lya emitting galaxies 
at high redshift have concentrated on applications related 
to probing the history of reionisation and the sources re- 
sponsible for that event. However recent attention has also 
focussed on obtaining large samples of Lya emitters to use as 



tracers of the density field (Hifl et al. 2004 20081, with appli- 
cation to cosmological distance measure and probes of dark 
energy at redshifts not previously studied (Hobby-Eberly 
Telescope Dark Energy Experiment, HETDEX). In this pa- 
per we discuss the implications of non-gravitational contri- 
butions to the observed clustering of Lya emitting galax- 
ies for studies of the shape, amplitude and angular depen- 
dence of the Lya galaxy power-spectrum. To this end we 
provide a simple framework that allows the influence of non- 
gravitational effects from Lya transmission to be evaluated 
with respect to the available precision of cosmological con- 
straints. As part of our study we determine the detail with 
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which the astrophysics governing the observed properties 
and flux of Lya emitters must be understood, in order for a 
large survey like HETDEX to achieve its theoretical perfor- 
mance in measurement of cosmological parameters. 

Our paper is set out as follows. In § [2] we construct a 
linear theory model for the power-spectrum of Lya selected 
galaxies, which is generalised to be applicable to any partic- 
ular model of Lya transmission. We next discuss a simple 
analytic model of transmission (§[3|, in which a fraction F of 
the Lya line is subject to optical depth tq. We provide an an- 
alytic formula for the resulting redshift space and spherically 
averaged power-spectra, which we use to discuss the physical 
origin of diflerent features in the power-spectrum. In §|4]we 
investigate the effect of transmission fluctuations on clus- 
tering using detailed, previously published models of Lya 
transmission. Having determined the likely extent of non- 
gravitational contributions to the observed clustering, we 
present an analysis of the likely impact of possible fluctua- 
tions in Lya transmission on the measurements of cosmolog- 
ical parameters in the HETDEX survey (§[5|. As a specific 



example we compute the Alcock-Paczynski effect ( Alcock & 
Paczynski 19791. In § |6] we turn this analysis around and 



discuss the precision that will be available for constraints 
on Lya transmission models. We present our conclusions in 
§ [7] In our numerical examples, we adopt the standard set of 
cosmological parameters ( Komatsu et al.|2011 1 , with values 
of fim = 0.24, fib = 0.04 and Qq = 0.76 for the matter, 
baryon, and dark energy fractional density respectively, and 
h = 0.73, for the dimensionless Hubble constant. 



2 MODEL FOR THE CLUSTERING OF LYa 
SELECTED GALAXIES 

We begin by briefly summarising the basic theory for ab- 
sorption of Lya photons from galaxies in the IGM, and then 
present a simple derivation of galaxy bias in linear theory. 
We use these as a basis to discuss the effects of fluctuations 
in the transmission of Lya flux through the IGM on the 
observed clustering of Lya emitters. Our derivation of fluc- 
tuations is similar to the analytic model presented in [Zheng] 
et al.|(|2011|). 



2.1 Lya absorption in the IGM 

The transmission of Lya photons (e.g. |Dijkstra et al.|2007ei| ) 
is summarised briefly here to provide context for the new 
clustering calculations. The total opacity seen by a photon 
initially at frequency 1/ is 



drnH{r)xH{r)aLy {1^ X [1 + V2{r)/c]) , (1) 



where the Lya absorption cross-section is written as a-Ly{v), 
and the gas at distance r greater than the virial radius r^ii 
has line-of-sight velocity Wz(r). In this expression xh is the 
fraction of hydrogen in atomic form, given at photoionisation 
equilibrium in the optically thin limit by 



nnar 



(2) 



where nn is the number density of hydrogen nuclii, F is the 
photoionisation rate and Orec is the Case-B recombination 



coefficient, a^ec = 4 .2 x 10-^^(Tgaa/10''K)-°-^cm^s-^ (e.g 
Hui fc Gnedin|1997[ ). Substituting we find 

n|(r)ai. 



1: 



dr 



r(r 



-rJLy {v X [l + Wz(r)/c]). 



(3) 



The total transmission is found by integrating over the fiux 
density J(i') in the Lya line 



r = 



(4) 



Approximating the Lya scattering cross-section as a 
delta function ai^y lv') ^ ai^y,tot5(v' — t^Ly), w here crLy.tot = 
JaLydu = (Rybicki & Lightman'l979' ^ and u'ir) = 

X [1 + V2{r)/c\, we can re- write equation (|3 1 as 



dr 



n^{r)arcc crLy.tot5(r — r^y) 

r(r) 



1 dv^ 
c dr 



(5) 



CnH(rLy)arccO-Ly 



r(rL: 



dv^ 



where v-^y — i/ll+Vz {r'Ly)/c] , in which rLy denotes the radius 
at which the photon is at resonance in the frame of the gas. 
Equation ^ illustrates the important point that the proba- 
bility of Lya scattering is proportional to the total number 
of hydrogens that the photon encounters per unit velocity 
along the line-of-sight. The optical depth is therefore pro- 
portional to the inverse of the line-of-sight velocity gradient 
(^^)~^, in addition to the more obvious dependencies of 
density squared and inverse of ionisation rate. Thus we find 



2 2-0.7(^-1) 

T OC 3 OC ^ . 

-prpO,^ dvx -p ' 

dr dr 



(6) 



where in the second proportionality we have assumed a 
power-law relation between density and temper ature of T oc 
with the polytropic index 7 = 1.4 ( |Hui fc Gnedin 



1997^ 



2.2 Galaxy bias and fluctuations in the number 
density of galaxies 

The likelihood of observing a galaxy at a random location 
is proportional to the local number density of galaxies. The 
likelihood of observing a galaxy within a region of large-scale 
overdensity is therefore equal to the ratio of the number den- 
sity of halos n{5) in a region of large-scale over-density 5 to 
the number density of halos in the background universe (n). 
This ratio has been used to derive galaxy bias for small val- 
ues of <5 ||Mo fc White|1996||Sheth et al.|2001| ). For example, 
in the Press & Schechter ( 1974 1 formalism we write 



dn{5)/dM 
dn/dM 



(i + S) 



l + <5 1 



dn 
dM 



+ 



d^n 


dv 




dn 


dMdv 


dS^ 




dM _ 



(j{M)v 



l + 5b. 



(7) 



where [dn/ dM)^^) and {dn/ dM){v) are the average and per- 
turbed mass functions, v = [5 — 1.69)/(t(M), a{M) is the 



^ Here, fa = 0.4167 denotes the oscillator strength for the Lya 
transition, and e and nie denote the charge and mass of the elec- 
tron, respectively. 



© 0000 RAS, MNRAS 000, 000-000 



4 Wyithe & Dijkstra 



variance in the density field smoothed on a mass-scale M at 
redshift z, and b is the bias factor. At small values of large 
scale overdensity 5, and in the absence of effects that influ- 
ence the relation between observed flux and halo mass, the 
number density of galaxies is proportional to [1 -I- 5b[M, z)]. 



2.3 The Lya emitter power spectrum 

In this section we estimate the effect of transmission fluc- 
tuations on clustering of Lya selected galaxies, beginning 
with equations Q and ([7| as motivation. Fluctuations in 
transmission of Lya radiation through the IGM modify the 
intrinsic luminosity that corresponds to an observed flux 
limit Fq. This modification of intrinsic luminosity in turn 
leads to a modification of number counts according to the 
luminosity function of Lya emitters. 

The density nLya of Lya emitters that are observed 
with fluxes greater than Fq [corresponding to an intrinsic 
luminosity Lq with transmission To = exp(— ro)] can then 
be expressed relative to the average nLyQ(> Lo, po, To) as 



WLyQ(> Fq) 

+(r- 



ni^ya{> Lo, po, To) X (1 + bS) 



To: 



dT 



dn^ya 



To.ro 



dT 



+{p- Po) 
+( 



or 

dT dni,ya 



Fo,To 



dT 
dT 



Fo.To 



a(arcom ) 



dn 



Lya 



To,Po 



dT 



(8) 



Fo,To 



where H is the Hubble parameter at scale factor a, and rcom 
is a co-moving distance. Here nLyQ(> Lo) is the mean num- 
ber density of Lya emitters with luminosities greater than 
Lo, and nLyc«(> ^^o) is the perturbed number density of Lya 
emitters with observed fluxes greater than the correspond- 
ing 7*0. We deflne the variable Sr = r/Fo — 1 as the fluctu- 
ation in the ionising background (discussed further below). 



We introduce a change of co-ordinates , 

and deflne the symbol 5^ = J^"' -j^ , which represents the 
fluctuation in line-of-sight velocity. With these we obtain 



n-Lyai> Fo) 



nLya(> Lo, Po, Fq) X (H- bS) 



+ 



9(logr) 



+ S 



91ogF 

Q(iogr) 

d log p 



dn 



Lya 



ro,ro a(logr) 

driLya 



Fo,To 



Fo,To 



a(iogr) 



9 log (dvz/dr) 



To ,po 



dULya 

d{logT) 



.(9) 



Fo.To 



Rearranging we get the fluctuation in the number- 
density of Lya emitters 



OLya 



^Lya 
fT^Lya 



1^5{b + Cp) + 5rCT + S^C^, (10) 



where 



Co 



1 dn{> L) 

flLya d{\ogT) 



7o,Lq 



djlogT) 
dlogp 



(11) 



To.Po 



and 



Cv 



dn{> L) 



rlLya 9(logT) 



a(iogr) 



5 log {dvz/dr) 



(13) 



To.ro 



By convolving the overdensity fleld of sources at points xo 
with a kernel oc exp [— (x — xo)/A]/(x — xo)^, the Fourier 
component of the ionising radiation fleld can be written 
( Morales fc Wyithe||2010[ ) 



5r{k) = bSk 



arctan (kX) 
kX ' 



(14) 



where A is the ionising photon mean-free-path. In redshift 
space the fluctuation in the density of Ly-a galaxies is 

dv. 



^Lya — ^Lya 



drcom Ha 



(15) 



Combining equations ([lOf, ( |14[ ) and ([15f, the Fourier com- 
ponent of the fluctuation in space density of Lya emitters is 
then 



l + Cr "'^°^7/^M +gp + (l^g^)/^' 



Slyaik) = 5k 

(16) 

wherfP^/ = d\og5/d\og{l + z), p is the cosine between the 
wave-number and the line-of-sight, and we have used the 
relation 5„(fc) — —fp?5k relating fluctuations in velocity 
gradient and density. The power-spectrum follows directly, 
yielding 

PLya{k,p)^ P{k) [b{l + CvKx)+Cp + {l-C^)fp^Y , (17) 

where P{k) is the mass power-spectrum and we have deflned 

arctan (fcA) 

(18) 



kX 



The spherically averaged power-spectrum is 



PZaik) 



= P{k)[{b{l + CrKx) + Cp)'' 

+ hbil + CrKx) + C,){l-a)f 



+ Hi-c.)'r 

5 



(19) 



We note that if Cr — Cp = Cv = (indicating no effect 
from Lya transmission) we obtain 



PLyaik,p)=PUk) [b + fp^] 



21 2 



and 



Plti'^) = P^ik) 



(20) 
(21) 



as expected in the standard case of galaxy clustering ( Kaiser 
|1987i ). 



3 ANALYTIC MODEL FOR THE EFFECT OF 
LYa TRANSMISSION ON CLUSTERING 

In this section we present a simple parameterised model for 
the fluctuations in Lya transmission, and use it to obtain 
analytic expressions for the constants Cp, Cr and Cv, and 



Cr 



dn{> L) 



riLya d{\ogT) 



d{\ogT) 



To,Lo 



aiogF 



To.ro 



(12) 



^ The quantity f is close to unity at high redshifts, taking values 
of 0.974, 0.988 and 0.997 at z = 2.5, 3.5 and 5.5. 
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hence for the power-spectrum. The resulting analytic ex- 
pression is instructive for elucidating the different important 
effects. In § [4] we use previously published detailed models 
of the Lya emitter transmission to calculate more physically 
motivated values for C„, Cr and Cv 



3.1 Simple model for Cp, Cr and Cv 

We begin by assuming a model in which the intrinsic Lya 
line of a galaxy is symmetric about the rest-frame Lya wave- 
length. We assume the Lya emitters to be located in an 
ionised IGM. We also assume that most absorption occurs 
in regions that are far enough from the galaxy that the ioni- 
sation rate is dominated by the ionising background, rather 
than by ionising flux associated with star-formation in the 
Lya emitting galaxy. In the absence of peculiar velocities in 
the IGM that cause departure from the Hubble flow, the red 
side of the line is transmitted through the IGM while the 



blue side is subject to resonant absorption (e.g. Madau 1995 
|Hu et al.'2004^ . However infall of intergalactic gas onto mas- 



sive galaxies leads to resonant absorption that can reach into 



the red side of the intrinsic line profile (e.g. Dijkstra et al, 

see 



2007a Laursen et al. 2011 



4.1 1. On the other hand, 



galactic outflows have the effect of redshifting the emergent 
Lya line relative to the true velocity of the galaxy, which 
makes the Lya photons more 'immune' to scattering in the 
IGM (see 



4.21 



To maintain generality we therefore assume that a frac- 
tion F of the line is subject to absorption in the IGM. As 
illustrated by equation Q, the transmission of Lya flux is 
dependent on the overdensity 5, and on the ionising back- 
ground and velocity gradient fluctuations 5r and Sv The 
remaining fraction (1 — F) is assumed to reach the observer. 
We therefore have the following expression for the overall 
transmission of the Lya line 



Ti5,5r,5v) = (1-F) + Fcxp 



-To- 



1 + ^(2.7-0.77) 



(22) 

where to is the mean Lya optical depth in the IGM as a 
whole. This model can describe a wide range of physical 
scenarios. The case in which the IGM suppresses only the 
blue half of the line corresponds to F — 0.5, while inf ailing 
intergalactic gas can result in F > 0.5. Conversely, scatter- 
ing through galactic winds can cause F < 0.5. 

To calculate Cp, Cr and Cv we first need an expression 
for dnrya{> L) / dlogT ■ For this calculation we assume that 
the luminosity function can be approximated as a power-law 
in the range of luminosities (L) observed. 



nrya{> L) — TlLya.O 



We then obtain 
dn{> L) 



dlogT 



dn{> L) 



dlosL 



L 



dlos L 



(23) 



r d log T 

Lq 

- l)nLyQ,0, 



where we have used the relation LT 
observed flux Fq). Utilising 



(24) 

I/o7o (i.e. constant 



dlogT 



dlogp 



7f),po 



-F(2.7 - 0.77)roe-^° 
(1 - F) -I- Fe-^o 



(25) 



and 

dlogT 
dlogT 

we obtain 



_ dlogT 

To.ro d log (dw^/dr) 



Froe-^" 
(1 - F) + Fe-^o ' 
(26) 



Cp^W- 1) 



(0.77 - 2.7)Froe^ 

[1- F) + Fe-^( 



(27) 



and 



Cr = C„ = (/3-l)- 



(28) 



Froe-^o 
' (1 - F) + Fe-^o • 

Putting these pieces together we write down an analytic 
estimate for the power-spectrum of Lya selected galaxies 



P„,(fc) [b + ffi'' + C [bKx + (0.77 - 2.7) - /m')] ' , (29) 
where 

(/3 - l)Froe-"« 



(1 - F) + Fe-^o ■ 
We also obtain the spherically averaged power-spectrum 
PlfLik) = P^{k)[{b + CibK, + {0.7^ ~2.7))f 



(30) 



+ -(6 + C(bA'A + (0.77-2.7)))(l-C)/ 



3.2 Results for analytic model 



(31) 



Some examples of the predicted power-spectrum of Lya 
selected galaxies calculated using this analytic model are 
shown in Figure [l] In each case the short-dashed, solid and 
long-dashed lines refer to the predicted power-spectrum as- 
suming the analytic model for Cp, Cr and Cv with F = 0.1, 
0.5 and 0.9 respectively. The power-spectrum in the absence 
of transmission effects (i.e. Cp = Cr = C„ = 0) is shown 
by the solid grey line. The upper sub-panels in each case 
show the ratio of these curves, indicating the fractional con- 
tribution of transmission effects to Lya emitter clustering. 
Three cases are presented in Figure [l] In the left column we 
show a source redshift of z = 5.5, and assume tq = 3 and 
A = lOOcMpc jBolton fc Haehnelt||2007| ), with a luminosity 
function slopa^ of /3 = 2. In the central column we show 
results for z = 3.0, tq = 1 and A = 300cMpc (Bolton & 



Haehnelt|2007| ), using the same luminosity function slope of 
/3 = 2. In the right hand column we again show results for 
z — 3.0, To = 1 and A = 300cMpc, but this time assume 
a steep luminosity function slope of /3 — 2.5. At each red- 
shift the simple analytic model predicts that fluctuations in 
transmission can lead to a factor of 2 difference or more in 
the power-spectrum amplitude. The effects are enhanced by 
a steep luminosity function, owing to the proportionality of 
the coefficients Cp, Cr and Cv to the factor (/3 — 1). 



^ The faint end slope of the Lyo luminosity function is not well 
constrained, and may shallower than ^ = 2 ( |Ouchi et ar]|2008| . 
However the LAEs that will be used in the HETDEX survey will 
be comparable in luminosity to those in the sample of |Ouchi et al.| 
l |2008[ l. At these luminosities, the slope of the luminosity function 
is steeper than at the faint end. 
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Figure 1. Clustering of Lya emitters in the simple analytic model. The Left, Central and Right panels show the the spherically averaged 
power-spectra in cases with {z,to,X,I3) = (5.5, 3, lOOMpc, 2), (^,ro,A, /3) = (3.0, 1, 300Mpc, 2) and (z,to,A,^) = (3.0, 1, 300Mpc, 2.5). In 
each case the short-dashed, solid and long-dashed lines refer to the predicted clustering assuming the analytic model for Cp, Cr and Cv 
with F = 0.1, 0.5 and 0.9 respectively. The clustering in the absence of transmission effects (i.e. Cp = Cr = Cv = 0) is shown by the 
solid grey line. Also shown (upper sub-panels) is the ratio of these curves, indicating the fractional contribution of transmission effects 
to Lya emitter clustering. The assumed halo mass was IO^^Mq. 



3.3 Contributions to the clustering amplitude 

We next calculate the difTerent contributions to the power- 
spectrum of Lya selected galaxies using the analytic model. 
The case of z = 3.0, ro = 1 and A = SOOcMpc, with lu- 
minosity function slope /3 = 2 is shown in Figure [2] The 
short-dashed, solid, and long-dashed lines refer to the pre- 
dicted power-spectrum assuming the analytic model in cases 
where contributions are included from Cp, from Cp and Cr, 
and from all of Cp, Cr and Cv respectively. We assumed 
F — 0.9 to accentuate the dependancies. The clustering in 
the absence of transmission effects (i.e. Cp = Cr = Cv ~ 0) 
is shown by the solid grey line. 

The simplest effect to understand is the decrease in the 
power-spectrum amplitude that follows the inclusion of den- 
sity fluctuations in the Lya transmission. The density fluc- 
tuations preferentially suppress observed flux in overdense 
regions, leading to lower observed galaxy number densities, 
and hence to a contribution to the clustering that counter- 
acts galaxy bias. Similarly, the inclusion of fluctuations in 
velocity gradient also decreases the power-spectrum ampli- 
tude. This is because fluctuations in the velocity gradient, 
which lead to increased transmission and hence to increased 
galaxy density at fixed Lya flux, are negatively correlated 
with over-density in mass. 

On large scales {kX <^ 1) the suppression of the power- 
spectrum by density fluctuations is counteracted by fluctu- 
ations in the ionising background. This is because the fluc- 
tuations in ionising background are both correlated with 
overdensity and biased, and so enhance transmission in 
overdense regions. In each example shown in Figure [ij the 
change of shape in the power-spectrum is evident at the scale 
k ~ 7r/(2A) corresponding to the ionising photon mean-free- 
path. On scales smaller than the mean-free-path {kX ^ 1), 
fluctuations in the ionising background are washed out and 
so do not contribute to modiflcation of the clustering am- 



F=0.9, z=3.0, To=1, A=300, /3=2 




Figure 2. Contributions to the spherically averaged clustering of 
Lya emitters in the simple analytic model. The case of z = 3.0, 
TO = 1 and A = 300cMpc, with /3 = 2 is shown. The short-dashed, 
solid and long-dashed lines refer to the predicted clustering as- 
suming the analytic model where contributions are included from 
Cp , from Cp and Cr , and fro all of Cp , Cr and Cv ■ We assumed 
F = 0.9 to accentuate the dependancies. The clustering in the 
absence of transmission effects (i.e. Cp = Cr = Cv = 0) is shown 
by the solid grey line. Also shown (upper sub-panels) is the ratio 
of these curves, indicating the fractional contribution of trans- 
mission effects to Lya emitter clustering. The assumed halo mass 
was 1O"M0. 



plitude. This is because only a fraction of the ionising back- 
ground is produced locally within the fluctuation, with the 
remainder being generated within a larger region that av- 
erages over many fluctuations. On these scales the fluctua- 
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Figure 3. Redshift space clustering of Lya emitters in the analytic model. The Central and Right panels show contours of the power- 
spectrum for the cases of {z, tq, A, j3) = (3.0, 1, SOOMpc, 2) and (z, tq, A, f)) = (3.0, 1, SOOMpc, 2.5). The Left panel shows the clustering in 
the absence of transmission effects (i.e. Cp = Cr = Cv = 0) for comparison. The solid contours are separated by a decade in Pj^y We 
have assumed F = 0.9, and a halo mass of IO^^Mq. 



tions in mass-density result in a lowering of the clustering 
amplitude by correlating lower than average Lya transmis- 
sion with overdensities of galaxies. It should be noted that 
our formulation ignores the possible Poisson contribution 
to fluctuations in the ionising background due to quasars. 
Poisson fluctuations introduce additional power beyond the 
component associated with the underlying density field of 
galaxies. The effect of Poisson fluctuations could become 
important at low redshift (z < 3), where quasars contribute 
significantly to the ionising background but is not included 
here. 



3.4 Redshift space clustering Lya galaxies 

In this sub-section we discuss the effect of the velocity 
structure in the IGM on the observed 2-dimensional red- 
shift space clustering of Lya selected galaxies on large 
scales. In Figure [3] we plot contours of the redshift space 
power-spectrum as a function of the line-of-sight (fcy = fc/i) 
and transverse {k±_ = k\/l — /x^) components of the wave- 
number k. The left panel shows the clustering in the absence 
of transmission effects (i.e. Cp = Cr = Cv = 0). The effect 
of infall in the linear regime ( Kaiser|1987 1 is clearly seen at 
small values of k, resulting in a power-spectrum amplitude 
that is increased at large scales in the line-of-sight direction. 

The central panel shows contours of the power-spectrum 
for the cases of {z, tq, A, /?) — (3.0, 1, 300Mpc, 2). We have as- 
sumed F = 0.9 and M = IO^Mq. The non-zero C„ term 



in this model counteracts the Kaiser (19871 effect, leading 



to more isotropic redshift space clustering on large scales. 
In the right panel we show a more extreme model for the 
luminosity function with (z, tq. A, f5) = (3.0, 1, 300Mpc, 2.5). 
In this case C„ > 1, and so the Kaiser (19871 effect is re- 



versed, leading to suppressed clustering transverse to the 
line-of-sight ( [Zheng et aL]|2011[ ). 



3.5 Coefficients in the Analytic Model 

The coefficients in the analytic model are dependent on as- 
sumed values for F and r. In Figure |4] we plot contours of 



the coefficients Cr and Cv (left panel), Cp (central panel) 
and the mean transmission To (right panel), each calculated 
as a function of F and tq. Large values of the coefficients 
require large values of F, indicating that a significant frac- 
tion of the intrinsic Lya line must be subject to absorption 
in order to influence the clustering of Lya selected galaxies 
at a level that is of order unity. 

It is interesting to ask what properties of the model are 
required to obtain clustering that is enhanced transverse to 
the line-of-sight as reported in the numerical simulations of 



Zheng et al. (20111. Inspection of equation (29 1 indicates 
that Cv > 1 is required. Assuming a luminosity function 
with P = 2, this can be achieved with F > 0.9 and to ~ 2 — 5, 
for which the transmission is To ~ 5 — 20%, in good agree- 



ment with the Lya emitters discussed in Zheng et al. ( 2011 1 



If the luminosity function is steeper, then the requirements 
on F are less stringent. For j3 = 2.5, the coefficients are 
increased by a factor of (/? — 1) = 1.5 so that C„ > 1 is 
obtained for F > 0.8. 



4 DETAILED MODELLING OF 

TRANSMISSION AND CLUSTERING OF 
LYa EMITTERS 

The analytic model discussed in §[3]is useful for investigating 
the qualitative dependencies of clustering in Lya selected 
galaxies. However a more detailed analyses is required to 
quantitatively predict the values of constants Cp, Cr and 
Cv, which describe the modification of the power-spectrum 
from that measured by a traditional galaxy redshift survey. 
In this section we describe calculation of Cp, Cr and Cv 
based on two previously published models of Lya emission. 
These models explore respectively, the effects of local star- 



formation and IGM infall (jPijkstra et al. 2007a I, and of 



galactic wind driven outflows (Verhamme et al.|2008 Dijk 



stra & Wyithe 20101, on the transmission of the Lya line 



through the circumgalactic IGM. 
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Figure 4. Contours of the coefficients Cr and C„ (left panel), Cp (central panel) and mean transmission To (right panel), calculated 
as a function of F and ro in our analytic model. We assumed (3 = 2 for these calculations. 



4.1 Modelling Lya transmission in the presence of 
IGM infall and galactic ionising flux 

The model presented in this section is based on the work in 



Dijkstra et al. (2007a I, to which the reader is referred for a 
full description of the calculations. The IGM transmission 
is calculated using a model for the IGM that accounts for 
clumping and infall. In this model, resonant absorption of 
LyQ photons by gas in the infall region (which extends out 
to several virial radii, see Barkana 2004) erases a significant 
fraction of the Lya line flux at frequencies redward of the 
LyQ resonance. Here, we briefly summarise the models main 
ingredients. 

The Lya flux from star-forming galaxies originates in 
the dense nebulae from which the stars form. Approximately 
two out of three (0.68) ionising photons produced by O stars 
(which are absorbed in the nebulae) are converted into Lya 
for case-B recombination (Osterbrock 19891. The total in- 



trinsic LyQ luminosity of a galaxy can then be estimated 
from 



I/Lya = 0.&9>hVa{l - /csc)Qh, 



(32) 



where hva = 10.2 eV is the energy of a LyQ photon, Qh 
is total luminosity of ionising photons, and /esc is the es- 
cape fraction of ionising photons from the galaxy. Since Qh 
depends on the number of O-stars, its value is sensitive to 
the assumed initial mass function (IMF) and metallicity of 
the gas from which the stars form. Under the assumption of 



constant star formation rate, Schaerer ( 2003 1 has calculated 



Qh for several different IMFs, and for a range of metallicities 
Z (expressed in solar units Zq). For a Salpeter IMF with 
lower and upper mass limits of M\ — IMq and Mu = IOOMq 
respectively 



logio(QH) = 53.8 + logio(MO - 0.0029(9 



-logio(^))'-', 
(33) 
"1 We as- 



in which M* is the star-formation rate in Mq yr 
sume Z = 0.05. For the models in this sub-section the shape 



4.2 



of the intrinsic Lya line is assumed to be Gaussian (see 
for a discussion of models where we assume different intrin- 
sic spectral line shapes). For gas that is optically thin to 
LyQ photons, a reasonable choice for the standard deviation 
of this Gaussian emission line is a a ~ Wvir, where Vvir is the 
virial velocity of the host galaxy halo ( Santos|2004 Dijkstra 



et al. 2007a I. We use equation (331 to calculate Qh as a 
function of . 

Given a galaxy spectrum blueward of the hydrogen ion- 
isation threshold {fa) of J(!^) oc u'^" , the photoionisation 
rate at distance r from the galaxy can be found from 



r(r) = Fbg + 



/3s O-oQufe: 



/3s 



47]- 7-2 



(34) 



where Fbg is the photoionisation rate of the meta-galactic 
background, and we have approximated the hydrogen pho- 
toionisation cross-section as anii^) ~ o-o{u/un)^^, with 
ao — 6.3 X 10~^*cm~'^. In addition to the radial depen- 
dence of the photoionisation rate we also require the radial 
density (p) and velocity (v) profiles of the intergalactic gas 
surrounding objects of mass M (,Barkana||2004[ ). Useful fit- 
ting formula ( Dijkstra et al.|2007a l for these are 



p(r) = 



20p(r/rv 
P 



r < lOrvir 
r > lOr-vir, 



(35) 



and 



l'infall(r) = 



H{z)r 



(r- 



where the velocity gradient 



Tvir < r < lOrvir 

r ^ lOrvir, 

(36) 

= (10rviri?(z) + 



'ycirc)/9rvir was choseu to make v{r) continuous. We eval- 
uate the optical depth in this model using equation (|4|, 
which is also integrated over the probability distribution of 
density contrasts ( [Miralda-Escude et al.||2000[ ). The radial 
dependencies of uh , v and F in equation (|4|) are specified by 
equations ( 34|36 1. 

For our calculations , we assume an I GM temperature 
of T = 2 X 10* K (e.g. |Lidz et al.||2010|), consider a halo 



( Orsi et al.|2008||Guaita et al.|2010[ ). 



mass of M = lO^'M© 

and a star-formation rate of M* = lOM0yr~^, which cor- 
responds to the mean and median UV derived SFRs for 
LAEs in the HETDEX Pilot survey (Blanc et al. 2011). 
The resulting total intrinsic luminosity in LyQ photons is 
2.3 X 10*=* [1 - /esc] erg/s. The escape fraction of 
ionizing photons, /esc, is uncertain (e.g. Yajima et al.|2010 
and references therein), and may vary significantly between 
individual objects (Shapley et al. 20061. We consider two 
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Figure 5. Example line profiles in the infall model, with modifications from adjustments in the ionising background F (left), the density 
p (center), and velocity gradient dv^/dr (right). The vertical dashed line indicates the Lya line center. Values for the sizes of the 
fluctuations considered in these quantities are listed in each case, together with the resulting transmission. The upper row shows results 
assuming that the galaxy does not contribute to the ionising flux. The lower row assumes an escape fraction of ionising photons from 
the galaxy of /esc = 10%. The assumed halo mass viras IO-'^-'^Mq. 



s ( Faucher-Giguere et al. 



cases in the papei[^ /esc = 0.0, and /esc = 0.1. We calculate 
profiles eX z = 3.0, for which Wcirc = 103 km/s and r-vir = 39 
kpc). The background photoionisat ion rate at z = 3.0 is 
taken to be Fbg = 0.5 x 10" 
[2008] ). 

Examples of the intrinsic (grey lines) and transmitted 
(black lines) Lya lines in this model are shown as a function 
of restframe wavelength in Figure [5] We show cases in which 
no ionising radiation escapes the galaxy (upper panels) , and 
in which the escape fraction is /oac = 0.1 (lower panels). The 
line profiles show the features of absorption redward of the 
intrinsic Lya central wavelength owing to the influence of 
gas infall, and no absorption redward of the blue most edge. 
The three sets of panels show the dependence of the line 
profile on fiuctuations in the density (p, left panels), ionis- 
ing background (Fbg, central panels), and velocity gradient 
(dv^jdr, right panels). For each quantity we present fluctu- 
ations (e.g. dp = p/po — 1, where po is the fiducial model) of 
±0.3 relative to the fiducial model. The resulting values of 
transmission T are presented in Figure [5] The transmission 
of the fiducial model is T = 0.42 for the case with /esc = 0, 
and T = 0.69 for the case with /esc = 0.1. 

We incorporate these fluctuations into our model as fol- 
lows. Firstly, to vary the density p, we make the modifica- 



The assumed /esc is degenerate with the assumed gas metallic- 
ity Z. 



tions nn nii{l + S), and (ii) Tgas T'gas(l + [7 — l]5) in 
equation (|3|. This temperature change affects the recombi- 
nation coefficient (see equatio n ^ . Second, for variations in 
r, we adjust Fbg in equation ( |34[ |. Finally, to study the im- 
pact of fiuctuations in the velocity gradient, we multiply the 
optical depth in equation ([T]) by a factor of (l-|-5„) in the lin- 
ear regime beyond lOrvir, and by a factor of (l-|-5infaii) in the 
infall region at r < lOrvir- Though not self-consistent, this 
procedure preserves the density and velocity profiles, and so 
isolates the effect of velocity gradient. To evaluate 5infaii we 
calculate fiuctuations in velocity gradient within the infall 

region relative to the fiducial model with . Noting 



that 



in the infall region, we keep the 



dr 9rvir 

circular velocity and virial radius fixed, but replace H with 
H {1 + 5v) to modify the velocity gradient beyond the infall 
region. This results in a fiuctuation in the velocity gradient 
within the infall region of 

_ dUinfall/d»" - dWinfall/d»"lo 



d^^infaii/drlo 

1 

1 -f- lOffr-vir/Uv 



(37) 



where in the last equality we have noted that the galaxy 
dynamical time is rvir/wvir ~ O.IH"^. Thus fluctuations in 
velocity gradient within the infall region are reduced relative 
to those in the linear regime. 

The line profiles in Figure [S] illustrate that fiuctuations 
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in density and ionising background lead to substantial modi- 
fication of the transmitted flux profile on the blue side of the 
Lya line, extending into the red side owing to infall. Fluctua- 
tions in density have a larger effect on the transmission than 
fluctuations in ionising background, and the effects have op- 
posite sign as discussed in § |3.3[ In the case where there 
is no contribution to the ionisation of the local IGM from 
the galaxy (/esc = 0), the effects of fluctuations in ionising 
background and velocity gradient are of similar magnitude 
in this model. However, in the case where galactic ionising 
flux also effects the ionisation state of the local IGM we flnd 
that the influence of fluctuations in ionising background is 
reduced. This is easy to understand since modiflcation of 
the ionising background level has little effect on the trans- 
mission in regions where the galaxy dominates the ionising 
flux. 

Based on these absorption proflles we can esti- 
mate the quantities flogT/c^logp, flogT/^logF and 
dlogT/dlog{dVz/dr), and hence the values of the constants 
Cp, Cr and Cv which govern the modiflcation of galaxy clus- 
tering. In the case where the galaxy makes no contribution to 
the ionising flux (/esc ~ 0), we flnd Cp = —0.72, Cr = 0.32, 
and Cv = 0.20. Inspection of Figure |4] shows that these val- 
ues are similar to those in our simple analytic model for 
F ~ 0.65 (corresponding to a fraction of the red half of the 
line having been absorbed due to infall), and r ~ 2 (which 
reproduces the infall model transmission at ~ 3). In the case 
where /esc ~ 0.1, the value of Cr = 0.05 is much smaller 
than for /esc = 0, owing to the reduced relative importance 
of the ionising background. Similarly, owing to the higher 
transmission, the /esc = 0.1 model also leads to lower values 



of Co 



-0.39 and Cv ~ 0.11. We can again flnd approxi- 



mate agreement when we compare these values of Cp and C 
to our simple analytic model (Figure |4| assuming F ~ 0.65 
and T ~ 0.5 (corresponding to the mean transmission of 
T = 0.7 for this model). 



Modelling Lya transmission in the presence 
of dust-free, symmetric ISM outflows 



4.2 



The modeling of § |4.1| assumes the intrinsic Lya line emerg- 
ing from the galaxy into the IGM to be a Doppler broadened 
Gaussian, which is symmetric in frequency around the Lya 
resonance. However, galactic outflows have the effect of red- 
shifting the emergent Lya line relative to the true velocity of 
the galaxy (e.g. Ahn et al.|2003 Verhamme et al.|2006 1, and 
there is strong observational evidence that this mechanism 
is at work. This evidence includes the observed blue-shift 
of interstellar metal absorption lines combined with the ob- 



served red-shift of the Lya emission line ( Steidel et al. 2010 1, 



column density A^hi, and outflow velocity Vsh- We assume 
that the shell has a radius of 1 kpc and a thickness of 0.1 kpc, 
but stress that the precise physical scale of the outflow is not 
important for our results. Our assumed gas temperature of 
TiSM ~ lO** K in the outflowing H I shell corresp onds to a b- 
parameter of ~ 13 km s^^ in the terminology of Verhamme 



et al. (20081. We further assume the H I shells to be dust- 



free (see § 4.3 for a discussion on dusty outflows). Verhamme 



et al.| (|2008|| typicaUy found that log Nm 

1 



19 - 21, and 
We therefore assume a model in which 
200 km s~^). We compute Lya 



Vah ~ — 500 km s 
(iVm,%h) = (10^" cm 
spectra emerging from the outflows using a Monte-Carlo 
transfer code ( Dijkstra et al.|2006 l. In our calculations, the 
Lya photons are emitted at line center (Arya = 1216A). We 
compute the impact of the IGM on the directly observed 
fraction of Lya by suppressing the intrinsic spectrum by 
exp(— r) (see § |4.1[ |. Further details on the calculation of 
this model can be found in [Dijkstra et al.| ( |2011[ ) . The grey 
solid line in Figure [5] shows an example of the Lya spec- 
tra emerging from the outflows. The emerging spectrum is 
highly asymmetric, with more flux coming out on the red 
side of the Lya line center. The spectrum peaks at about 
~ 2wsh, as expected for radiation that scatters back to the 
observer on the far side of the galaxy (see |Ahn et aL]|2003| 
[Verhamme et aL]|2006| for a detailed discussion on these 
features in the spectrum). 

In Figure |6] we show the transmitted (black lines) Lya 
line for the outflow model. As in Figure [5] we show cases in 
which no ionising radiation escapes the galaxy (upper pan- 
els), and in which the escape fraction is /esc = 0.1 (lower 
panels). As before we show the dependence of the line pro- 
flle on each of fluctuations in the density (p left panels), 
ionising background (Fbg, central panels), and velocity gra- 
dient {dvz/dr, right panels). The resulting values of trans- 
mission T are listed. The transmission of the flducial out- 
flow model is T = 0.85 where /esc = 0, and T = 0.94 where 
/esc = 0.1. These values are larger than were found for the in- 
fall case, owing to the large fraction of radiation that scatters 
away from resonance before emerging from the galaxy. All 
trends of the transmission with fluctuations in density, ion- 
ising background and velocity gradient previously described 
for the infall model, are also present in the outflow mod- 
els. However the modifications of the transmitted flux are 
smaller than are found in the infall case. This is because a 
smaller fraction of the emergent flux is subject to resonant 
absorption in the IGM. 

Based on these absorption proflles we can estimate the 
values of the constants Cp, Cr and Cv which govern the 
modiflcation of galaxy clustering in the outflow model. In 
the case where the galaxy makes no contribution to the ion- 



and the fact that Lya line shapes are asymmetric at all red- 
shifts (e.g. [Mas- Hesse et al.|2003||Heckman et al.|2011[ ). We 
refer the reader to Dijkstra et al. ( 2011 1 for a more extended 



discussion. Verhamme et al. ( 2006 2008 1 have developed a 



ising flux, we flnd Cp 



-0.06, Cr = 0.04 and C„ = 0.02. 



Inspection of Figure [4] shows these values are similar to our 
simple analytic model where F ~ 0.2 (corresponding to most 
of the line having been scattered redward of Lya by the out- 
flow), and r ~ 2 (which results in a transmission equal to 
that predicted by the detailed outflow model). In the case 
where /esc = 0.1, the value of Cr = 0.005 is much smaller 
than for the /esc = case owing to the reduced importance 
of the ionising background. Similarly, because of the higher 
transmission, the model with /esc = 0.1 also leads to lower 
and Dijkstra & Wyithe (20101, we model the outflow as a values of Cp = —0.05 and Cv = 0.01. We again flnd approx- 
spherically symmetric thin shell of gas that contains an H I imate agreement when we compare these values for Cp and 



simple model in which scattering of Lya photons by H I in 
these outflows successfully explaines the observed Lya line 
shapes observed in Lya emitting galaxies at 2: = 3 — 6 (also 
see Vanzella et al.|2010 |. 

In this section we repeat the exercise of §|4.1|for a suite 



of outflow models. Following Verhamme et al. ( 2006 2008 1 
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Outflow model /pop = 



T 1 1 1 — I I — I 1 1 1 — I I — I 1 1 r 




Outflow model /per = 0.1 




Figure 6. Example line profiles in the outflow model, with modifications from adjustments in the ionising background F (left), the 
density p (center), and velocity gradient dv^/dr (right). The vertical dashed line indicates the Ly« line center. Values for the sizes of the 
fluctuations considered in these quantities are listed in each case, together with the resulting transmission. The upper row shows results 
assuming that the galaxy does not contribute to the ionising flux. The lower row assumes an escape fraction of ionising photons from 
the galaxy of f^sc = 10%. The assumed halo mass was IQ-'^-'^Mq. 



Cv to our simple analytic model (Figure |4]| with F ~ 0.2 
and r ~ 0.5, corresponding to the mean transmission of 
T = 0.94 for this model. 



4.3 Modeling the IGM transmission for dusty, 
anisotropic outflows 



The calculations presented in § |4.2| ignore the impact of dust 
on the Lya radiation field. However dust can play an im- 
portant role in the scattering of Lya photons within galax- 
ies (e. g. |Neufeld|[l99l| [Hansen fc Oii||2006| . [Laursen etldT 
( |2009[ ) have performed Lya radiative transfer calculations 
in simulated galaxies, finding that the effect of dust is to 
narrow the Lya line emerging from a galaxy relative to the 
dust-free casa^ Narrowing the Lya line causes a larger frac- 
tion of Lya photons to emerge at frequencies where they 
are subject to scattering in the IGM. We therefore expect 
the impact of the IGM to be stronger in cases where dust is 
included in the modelling of outflows. 



^ This is mostly because Lyo scattering can be described by a 
diffusion process in both real and frequency space (see |Dijkstra| 
|et al.|[2006[ and references therein). A large displacement in fre- 
quency requires a large number of scatterings, and therefore a 
long trajectory through the scattering medium. As the dust con- 
tent of this medium is increased, the probability that the photon 
is destroyed by a dust grain is enhanced. 



In our model, the galactic outflow is represented by a 
spherical shell. Departures from this idealized gas distribu- 
tion should also lead to a larger impact of the IGM on the 
emerging Lya line. This is because Lya photons will escape 
more easily from outflows in which either the covering factor 
is less than unity (because some sight-lines simply do not in- 
tersect with the outflowing material) , or when the scattering 
medium is clumpy. The latter is demonstrated by |Hansen fc| 
Oh ( 2006 1 who have studied Lya transfer through clumpy 



outflows, and have shown that a fraction of the photons can 
escape at line center. These results suggest that more com- 
plicated, and realistic models of winds than those employed 
in this paper will result in a stronger impact of the IGM on 



the observed Lya flux than we have computed here (§ 4.2 
also see Barnes et al.|20lT I . 



4.4 Summary of detailed modelling 

The values of Cp, Cr and Cv found from the detailed mod- 
elling in this section are summarised in Table [l] The pre- 
diction from the infall model for Lya transmission is that 
that there will be signiflcant contributions to the observed 
power-spectrum from fluctuations in the Lya transmission. 
Indeed, if the escape fraction of ionising photons is very 
small, then the contributions are expected to be of order 
unity. In contrast, the prediction from the outflow model is 
that contributions to the observed power-spectrum will be 



© 0000 RAS, MNRAS 000, 000-000 



12 Wyithe & Dijkstra 



Table 1. Evaluations of the constants Cp, Cr and Cv for the 
different models. 



model 




Cp 


Cr 




infall (/osc = 


= 0) 


-0.72 


0.32 


0.20 


infall (/osc = 


0.1) 


-0.39 


0.05 


0.11 


outflow (/esc 


= 0) 


-0.06 


0.04 


0.02 


outflow (/esc = 


= 0.1) 


-0.05 


0.005 


0.01 



an order of magnitude smaller, at the level of ~ 5 — 10%, 
although these numbers are likely to be conservatively small 
(see 



4.31 



These results suggest that measurement of the terms 
Cp, Cr and Cv from the observed power-spectrum of Lya 
selected galaxies will provide a new avenue to study the 
relationship between the Lya flux of galaxies and their local 
IGM. On the other hand, the expected non-zero values of 
Cp, Cr and Cv may complicate attempts to use the power- 
spectrum of Lya selected galaxies to constrain cosmological 
parameters. We turn to this topic for the remainder of the 
paper, in which we present an application of our general 
model for clustering of Lya selected galaxies to the planed 
HETDEX survey. 



5 LYa TRANSMISSION FLUCTUATIONS IN 
GALAXY REDSHIFT SURVEYS 

One of the primary science drivers motivating large galaxy 
redshift surveys is measurement of dark energy, and its evo- 
lution. Traditional galaxy redshift surveys are best suited to 
studies of the dark energy equation of state at relatively late 
times {z < 1) due to the difficulty of obtaining accurate red- 
shifts for a sufficiently large number of high redshift galaxies. 
Although detection of the Integrated Sachs- Wolfe effect puts 
some constraints on the integrated role of dark energy above 
z ~ 1.5 [see, e.g., [Giannantonio et al.| ( |2008[ ) and references 
therein], we currently have very limited information about 
the nature of dark energy at high redshift. If dark energy 
behaves like a cosmological constant, then its effect on the 
Hubble expansion is only significant at z < 1 and becomes 
small at z > 2. In this case, studies of the power-spectrum at 
low redshift would provide the strongest constraints. How- 
ever, as the origin of dark energy is not understood we can- 
not presume a priori which redshift range should be studied 
in order to provide optimal constraints on proposed models. 
Probes of dark energy at higher redshifts have been sug- 
gested. These include measurements of the power-spectrum 
from a Lya forest survey, which could potentially be used 
probe the evolution of dark energy through measurement 
of the baryonic acoustic oscillation (BAO) scale scale for 
redshifts as high as z ~ 4 (McDonald &: Eisenstein 2007| ). 
Similarly, studying the temporal variation of high resolution 
quasar spectra may probe the evolution of dark energy in the 
window 2 < z < 5 ( [Corasaniti et al.|2007[ ). 

The Hobby-Eberly Telescope Dark Energy Experiment 



( Hill et al.||2004[ |2008[ HETDEX), promises to provide a 
very important advance in our understanding of dark en- 
ergy, by measuring its contribution to the energy density 
at high redshift {z ~ 2.5) where there are currently no di- 



rect constraints. At the same time, a precision measurement 
of curvature will assist in breaking the degeneracy between 
dark energy and curvature present in lower redshift experi- 
ments ( |Hill et al.|[2004} [McDonald fc Eise nstein|[2007| . The 
HETDEX approach is to obtain approximately 0.8 million 
redshifts of Lya selected galaxies at 1.9 < z < 3.5. These 
galaxies will be obtained over an area of 400 square de- 
grees, with a survey volume of V ~ OGpc"^, and a galaxy 
space density of n^y ~ 10~*Mpc~"^. In the absence of non- 
gravitational contributions to the clustering, such a survey 
is able to use the measured power-spectrum to determine 
the local Hubble expansion at z ~ 2.5, and the angular di- 
ameter distance out to z ~ 2.5 to 0.8% each. In this section 
we discuss the influence that the non-gravitational contribu- 
tion to the observed power-spectrum of Lya selected galax- 
ies may have on the precision of cosmological constraints 
from a survey like HETDEX. We find that the transmission 
can strongly influence the constraints that are available, and 
determine the precision with which the parameters Cv, Cr 
and Cp will need to be understoc 
achieve its theoretical precision. 



5.1 Power-spectrum and power-spectrum 
sensitivity 

In this section we describe the power-spectrum, and the esti- 
mate of power-spectrum sensitivity, that we employ to calcu- 
late the influence of transmission fluctuations on cosmologi- 
cal constraints that will be available in HETDEX. Analysis 
of the galaxy power-spectrum derived from N-body simu- 
lations has shown ( Seo fc Eisensteiri||2005 1 that the power- 
spectrum can be treated as linear on scales greater than 
15 co-moving Mpc (i.e. fcmax ^0.4Mpc~^) at z = 3.5, in- 
creasing towards higher redshifts. Our assumption of a linear 
mass power-spectrum should therefore be sufficient for this 
analysis. However, weak oscillatory features in the power- 
spectrum, such as the baryonic acoustic oscillations, are sup- 
pressed on even larger scales because matter moves across 
distances on the order of ~ 5-lOMpc over a Hubble timfQ 
As groups of galaxies form, the linear-theory prediction for 
the location of each galaxy becomes uncertain, and as a 
result noise is added to the correlation among galaxies and 
hence to the measurement of the mass power-spectrum. The 
noise associated with the movement of galaxies smears out 
the acoustic peak in the correlation function of galaxies in 
both real and space ( [Eisenstein et al.|2007| [Seo et al.|2008| ). 
The associated reduction of power in the baryonic acoustic 
oscillations is found to be in excess of 70% on scales smaller 
than fcniax ~ 0.4Mpc~'^ at z ~ 3, corresponding to a length 
scale of ~ 7r/(2fcniax) = 3.9 comoving Mpc ( Seo et al.|2008 |. 
For our cosmological analysis we therefore use the following 
mass power-spectrum 



(fc) 

+ (P„,(fc) -P,n,nw(fc))exp I -fc^^^ ^ 



,(38) 



^ This characteristic scale of displacement follows from the fact 
that CTg, the normalisation of the power-spectrum on 8/i^^Mpc, 
is of order unity at the present time. 
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HETDEX z=2.5 V=9 Gpc^ n„a|=10"''Mpc"3 £i+C. = 2.5 



10' 
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k (hMpc^) 
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k (h Mpc"'') 



0.3 



Figure 7. Left panel: The predicted spherically averaged power-spectrum of Lya galaxies, using equation | |19| l with the combinations 
[Cv, bCr, {b + Cp)] = (0, 0, 2.5) and (0.5, 0.1, 2.5). Right panel: The baryonic acoustic oscillations, plotted as fluctuations relative to the 
"no wiggle" power-spectrum. The error-bars assuming the observational survey parameters V ~ 9Gpc'^ and nLy ~ 10~^Mpc~^ that are 
appropriate for HETDEX, with each point providing an independent measures of the power. 



where Pm.nw is the "no wiggle" form from lEisenstein fc| 5.2 



Hu| |T999| ). The non-linear scales in this expression are 
Sx = 4.6[(l-H2)/3.5]-^ Mpc and Sy = 9.2[(l-f zj/S.Sj-^Mpc 
( Seo fc Eisenstein|2007| in the high redshift limit. Note that 
we consider only scales k < 0.4Mpc~^, and so do not include 
the effects of the fingers-of-god that arise from random mo- 
tions within virialized halos in our analysis. As shown in 



Shoji et al. (2009 1, this has no influence on the cosmological 



constraints inferred from the large scale power-spectrum. 

Like traditional galaxy redshift surveys, the observed 
LyQ galaxy power-spectrum is sensitive to the underlying 
mass power-spectrum (Pm) and galaxy bias 6. However in 
addition, there is also dependence on the parameters Cr, 
Cp, Cv, which are related to properties of Lya transmission 
through the IGM, and on the mean-free-path A. Inspection 



of equations (17 1 and (19 1 indicates that not all of 6, Cp 



and Cr can be measured independently. Instead, the power- 
spectrum of LyQ emitting galaxies depends on the param- 
eters bCr, (Cp + b) and Cv In the left panel of Figure [7] 
we plot the predicted spherically averaged power-spectrum 



of Lya galaxies, using equation ( 19 1 with the combinations 
\Cv,hCT,{h + Cp)] = (0,0,2.5) and (0.5,0.1,2.5). These pa- 



rameters are motivated by an extreme outflow model (§ 4.2 1, 
and an infall model ( 



4.1 1 respectively. 



Figure [7] also presents error-bars assuming the observa- 
tional survey parameters V ~ 9Gpc^ and TiLy ~ 10~'*Mpc~'^ 
appropriate for HETDEX. We evaluate the uncertainty in a 
k-space volume 27rfe^ AA: sin(/i) A/i as 



APLy = Pi 



Ly 



PLynLy 



k Afc sin(/i)A/i 



V 



(2^) 



-1/2 



(39) 

where the sum in the first term encapsulates cosmic variance 
and galajcy shot-noise respectively, and the second term cor- 
responds to the number of modes measured in the survey. 
We find that shot-noise dominates at fc > O.lMpc"^. The 
right panel of Figure [T] shows the baryonic acoustic oscil- 
lations, plotted as fluctuations relative to the "no wiggle" 
power-spectrum. 



Measurement baryonic acoustic oscillations 
with Lya selected galaxies 



The imprint of BAOs on the mass power spectrum provides 
a cosmic yardstick that can be used to measure the depen- 
dence of both the angular diameter distance and Hubble 
parameter on redshift. The wavelength of a BAO is related 
to the size of the sound horizon at recombination. Its value 
depends on the Hubble constant, and on the dark matter 
and baryon densities. However, it does not depend on the 
amount or nature of dark energy. Thus measurements of 
the angular diameter distance and Hubble parameter can in 
turn be used to constrain the possible evolution of dark en- 
ergy with cosmic time (e.g. Eisenstein et al.|1998 Eisenstein 
|2002). 

Importantly, a measurement of the BAO scale is not 
subject to modifications of the overall shape or angular de- 
pendence of the power-spectrum, owing to non-linear grav- 
itational growth. To illustrate this point |Seo fc Eisenstein] 
(2005) modelled the power-spectra from a series of N-body 
simulations using the addition of a linear power-spectrum 
and a scale dependent polynomial to describe galaxy bias 



and anomalous power. Seo & Eisenstein ( 2005 1 find that 



they are able to recover the BAO signal by subtracting a 
smooth function from the matter power-spectrum measured 
in their N-body simulations. 

Our model is intrinsically linear and so does not include 
scale dependent bias or anomalous power. However the sub- 
traction of a smooth function to recover the BAO signal 
is valid for a range of different scale dependent contribu- 
tions (e.g. |Mao et al.|2008|[Rhook et aL]|2009[ ). As a resuh, 
while Lya transmission fiuctuations inhibit the extraction 
of the full cosmological information through consideration 
of the whole power-spectrum, they should not significantly 
alter the ability of a Lya galaxy survey to measure the BAO 
scale. 

To illustrate this point we have fitted the analytic ap- 
proximation to the baryonic oscillation component of the 
redshift space power spectrum following |Glazebrook fc Blake| 
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Figure 8. Constraints on the line-of-sight and parallel BAO scale. The three panels show contours of likelihood for the parameter set 
{a±, ct\\), in the cases of a traditional galaxy power-spectrum (left), and Lya galaxy power-spectra with [Cv, bCr, {b + Cp)] = (0, 0, 2.5) 
(central) and [C„, bCr, {b+Cp)] = (0.5, 0.1, 2.5) (right) respectively. Contours of likelihood for the values of ay and a± are shown (at 61% 

and 14% of the peak likelihood). The values of 61% ~ e~^^ 1'^ and 14% ~ e~^^ are chosen to equal the contour height corresponding 
to Ax^ = 1 and Ax^ = 4, and therefore the projection of these contours onto the axis for a particular parameter represents the 1-sigma 
and 2-sigma ranges respectively. The 61% contour for the traditional galaxy redshift survey constraints is repeated in the central and 
right panels for comparison (thick dotted line). 



(20051 



P\,y(k\\,k^_) = PLy,nw(fc||,fc±) 

X I + Ak cxp ' ' ^ 



O.OTMpc"^ 



X sin 2tv 



{l + aj_)kA 
X exp I —k 



(1 + a||)fcA, 



2Si(l-M') + Si|M' 



,(40) 



to estimate the constraints on the hne-of-sight and trans- 
verse BAO scales (ay and a±). In this expression the "wig- 
gle free" power-spectrum (PLy,nw) is computed using equa- 
tion (171, with the mass power-spectrum (Pm) replaced by 



the "wiggle- free" mass-power-spectrum (Pm.nw)- The ob- 
served power-spectrum Pi^y is modelled as the sum of PLy.nw 
and a decaying sinusoid with characteristic periods in the 
line-of-sight and transverse directions of (1 -|- a||)fcA and 
(1 + a±)kA- We include the factor of non-linear suppres- 



sion of the BAO amplitude (Seo & Eisenstein' 2007). This 



function has three parameters A, ay and a±. The value of 
A is determined to high accuracy from observations of the 
Cosmic Microwave Background. For the purposes of this 
analysis we therefore assume that A is a known constant 
(namely A — 2.1), and fit only for Qf|| and a± (around the 
best fit value of in the absence of noise). We fit only to 
values of fc < 0.4Mpc~^. With this parameterisation, the ac- 
curacy with which ay and q_l can be measured determines 
the constraints that BAO can place on the line-of-sight and 
transverse distances, and hence on the Hubble parameter H 
and angular diameter distance Da respectively. 

In Figure|8]we show contours of a|| and a± derived from 
the BAO analysis. For graphical representation of our results 
we show contours at 61% and 14% of the peak height. The 
values of 61% ~ e"'^''^^^ and 14% ~ e'^^^^'^ are chosen 
to equal the contour height corresponding to Ax^ = 1 and 
Ax^ = 4, and therefore the projection of these contours onto 



the axis for a particular parameter represents the 1-sigma 
and 2-sigma ranges respectively. The left hand panel shows 
the expectations for a galaxy redshift survey with parame- 
ters corresponding to HETDEX. Our analysis yields 1-sigma 
errors of Aan ~ 3.0% and Aa± ~ 2.0%. For consistency we 
compare with the BAO constraints for HETDEX presented 
in Shoji et al. (20091. These authors find similar constraints 
of Aay ~ 2.5% and Aqx ~ 1.8%. 

In the central and right hand panels of Figure [8] we 
allow for scale and direction dependent modifications to 
the power-spectrum owing to fluctuations in Lya transmis- 
sion (equation |17[ ). In these cases the fiducial models have 
(C„,feCr) = (0,0) and (C„,&Cr) = (0.5,0.1) respectively. 
No prior probabilities on their values were assumed. There 
is a very small difference in the resulting BAO constraints, 
indicating that modification of the power-spectrum by fiuc- 
tuations in Lya transmission will not inhibit use of the BAO 
scale for studies of dark energy in Lya selected galaxy sur- 
veys. The values for constraints on ay and a± are listed in 
Tabled 



5.3 Application of the Alcock-Paczynski test 



Shoji et al. ( 2009 \ have argued that much more accurate 



constraints on cosmological distances are obtained by con- 
sideration of the whole power-spectrum shape rather than 
only the BAO scale. To quantify the potential of the Lya 
galaxy power-spectrum for measuring cosmological parame- 
ters, we therefore calculate the Alcock-Paczynski effect ( Al- 
cock fc Paczynski|1979[ ). Our approach is to specify the gen- 
eral result of Barkanay ( |2006[ ) for the distortion of the true 



power-spectrum [P^y{k, ^)] that results from an incorrect 
choice of cosmology. Dilation parameters a and a± are used 
describe the distortions between the transverse and line-of- 
sight scales, and in the overall scale, respectively. These are 
defined such that (1 -I- a) is the ratio between the assumed 
and true values of (DaH), while (1 -f- a±) is the ratio be- 
tween the assumed and true values of the angular diameter 
distance, Da- In the Alcock-Paczynski test, the correct cos- 
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Figure 9. Constraints on power-spectrum distortions achievable via the Alcock-Paczynski test. The three panels show contours of 
likelihood for the parameter set {a±, <^\\), in the cases of a traditional galaxy power-spectrum (left), and Lya galaxy power-spectra with 
[Cv,bCT,ib + Cp)] = (0,0,2.5) (central) and [Cv , bCr , (b + C p)] = (0.5,0.1,2.5) (right) respectively. Contours of likelihood for the values 
of aj_ and Q|| are shown (at 61% and 14% of the peak likelihood). The 61% contour for the traditional galaxy redshift survey constraints 
is repeated in the central and right panels for comparison (thick dotted line). 



mology is inferred by finding cosmological parameters for 
which a = a± = 0. 

To calculate the Alcock-Paczynski effect we apply equa- 



tion (8) in |Barkana| ( |2006 l 

PLy(fc,At) = (l + a — Saj 



+a(l- 



(41) 



to the power-spectrum in equation (171. Here the power 



spectra Phy and Pi^y are evaluated at the observed k. This 
procedure results in a modified power-spectrum that is re- 
lated to the true mass density power-spectrum (P^) via 

PLy(fc,M) = PL{{b + Cp) + bCrK{k) + {l^C^)fi^f 
(1 + a - 3q±) + -77-7^ (qm^ - QJ 



+ 



dlnfc 

PL ((6 + Cp) + bCrKik) + (1 - 



n dK 2 



- Qx) 



+4a(l - a) (1 - /i") /m'] • (42) 

For dark energy studies it is more interesting to con- 
strain the quantities Da and H independently rather than 
Da and the product DaH. We define the dilation parame- 
ter an such that (1 -|- ay) is the ratio between the assumed 



values of H. The dilation parameter a in equation ( 42 1 is 
then expressed as 



1 + a|| 



(43) 



Equation ( 42 1 can then be used to find the precision of con- 



straints on H and Da- Inspection of equation (421 suggests 



that some degeneracies are expected between cosmologi- 
cal constraints (parameterised by ay and ax) and the un- 
known parameters describing the modification to the power- 
spectrum owing to Lya transmission fiuctuations. 



Alcock-Paczynski constraints on the 
power-spectrum 



5.4 



We next use equation (|42| to calculate the permissible 
region of parameter space p = (a||,a_L,6 -|- Cp,C„,feCr) 
around a true solution with power-spectrum P^y and po = 
(0,0,2.5,0,0). We have assumed that the mean-free-path 
of ionising photons is known a'priori, and do not fit it as 
a free parameter. We take the value to be Amfp = 300 
co-moving Mpc ( [Bolton fc Ha ehnclt 2007', Faucher-Giguere] 
et al.|2008[ ). Using the power-spectrum sensitivity specified 
in equation (39 1, we construct likelihoods 



\nC(p) = 



2^ 



APLy(fc,M) 



-f In Cc„ + hi £c„ + In Cc^ 



(44) 



where the sum is over bins of k and /i, and Ccp, Cc^ and 
Ccj. are the a-priori likelihoods for the parameters b + Cp, 
Cv and hCr respectively. To account for the possibility of 
non-linearity in the smooth power-spectrum at small scales 
we restrict our fitting to wave numbers fcmax < 0.4Mpc^^. 



5.5 Constraints for a traditional galaxy redshift 
survey 

To provide a baseline for our analysis we first consider 
the case where Lya transmission has no effect on observed 
galaxy flux (i.e. we consider only power-spectra with Cp = 
Cv ~ Cr = 0), as is the case for a traditional galaxy redshift 
survey. For this analysis we further assume that the value 
of / is known and that the shape of the primordial power- 
spectrum is well measured by other means. We refer to these 
constraints as being for a traditional galaxy redshift survey 
in the remainder of this paper. As discussed in § |5.8[ |Shoji| 



et al. ( 2009 1 have investigated the consequences of relaxing 



these strict prior constraints 
As noted in |5.3 



such that (1 + a|i7 
true values of (//), while (1 + a 



the parameters ay and a± are defined 
is the ratio between the assumed and 
is the ratio between the 



assumed and true values of the angular diameter distance, 
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Figure 10. Constraints on power-spectrum distortions and transmission models acliievable via the Alcock-Paczynski test. The six panels 
show contours of likelihood for the parameter sets (Ct,,a||), {aj_,Cv), {a±,bCY), {6Cr,«||), {C'v,b + Cp) and {bC'r , b + C p) , in the case 
of a hya galaxy power-spectrum with [Cv, bCr, {b + Cp)] = (0, 0, 2.5). Contours of likelihood for the values of ay and a± are shown (at 
61% and 14% of the peak likehhood). 



Da- Thus, the precision with which a± can be measured 
provides an estimate of the relative precision with which Da 
can be measured. Similarly, precision with which the local 
value of H can be measured is provided by the precision 
with which aii can be measured. 



The left hand panel of Figure [9] presents likelihood con- 
tours for the parameter set (Q:^,a||) obtained from a tra- 
ditional galaxy redshift survey assuming the HETDEX vol- 
ume and galaxy density. The likelihood is marginalised over 
the bias 6 assuming a flat prior probability. The contours 
show a degeneracy between an and qj_. This negative corre- 
lation is fundamental to the Alcock-Paczynski effect. When 
the redshift space distortion is known perfectly well, the 
departure of the real-space power spectrum from isotropy 



can be used to determine DaH (or a in equation 41 1. Fig- 
ure [9] indicates that line-of-sight and angular distortions of 
the power-spectrum compared with an assumed model can 
be measured at better than the ~ 1% level, indicating that 
our analysis is consistent with expectations for the HET- 
DEX survey. Comparison with Figure |8] shows that the pre- 
cision available on the line-of-sight and radial distances mea- 
sured from the Alcock-Paczynski test using the full power- 
spectrum shape are a factor of several better than from an 
analysis of the BAO scale alone ( Shoji et al.||2009 l. Values 
for these and subsequent constraints on a|| and ax are listed 
in Table [21 



5.6 Constraints for a Lya galaxy redshift survey 

We next calculate the level to which Lya transmission fluc- 
tuations influence the precision with which ay and a± can be 
measured. The central panel of Figure [9] presents likelihood 
contours for the parameter set (aj_ , Cfy ) obtained from a Lya 
galaxy survey, again assuming the HETDEX volume and 
galaxy density, with a fiducial model having [Cv,bCr, {b + 
Cp)] = (0,0,2.5). Here the likelihood is marginalised over 
the parameters Cv, bCr and (6 -I- Cp) assuming flat prior 
probabilities (i.e. Ccp ~ 'Cc„ = jCc^ = const). The 1-sigma 
contour for the traditional galaxy redshift survey case is 
repeated for comparison. Figure |9] indicates that without 
knowledge of the detailed properties of Lya transmission 
(enabling prediction of Cv, Cr and Cp), the line-of-sight and 
angular distortions of the power-spectrum can be measured 
at the ~ 1.3% level, a factor of ~ 1.5 decrease in the avail- 
able cosmological precision relative to a traditional galaxy 
redshift survey. 

To calculate the constraints on ay and a± for a Lya 
galaxy survey we must specify the transmission model. In 
the above calculation we assumed a fiducial model with no 
transmission effects, but allowed for their possible existence 
when performing the cosmological fit. In the right panel of 
Figure [9] we repeat this analysis for an assumed fiducial 
model which has [Cv,bCr, {b + Cp)] = (0.5,0.1,2.5), and so 
includes strong modification of the observed power-spectrum 
from fiuctuations in Lya transmission. In this case we find 
that constraints on the angular distortions of the power- 
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Figure 11. Constraints on power-spectrum distortions and transmission models achievable via the Alcock-Paczynski test. As per 
Figure [lO| but for the case of a Lya galaxy power-spectrum with [Cv, feCr, {b + Cp)] = (0.5, 0.1, 2.5). 



spectrum are unchanged (~ 1.3%), but that the line-of-sight 
distortions can only be measured at the ~ 1.7% level. 

To investigate the origin of the decrease in precision 
that follows inclusion of Lya transmission fluctuations in 
a power-spectrum analysis, we show contours of likelihood 
for the parameter sets (Ci,,Q||), (ax,C„), (a±,bCr) and 
(6Cr,a||), in the cases of [Cv,bCr, {b + Cp)] = (0.5,0.1,2.5) 
and [C„,6Cr,(6 + Cp)] = (0,0,2.5) (Figures [lO] and 



11 



spectively). For each set the likelihood is marginalised over 
the remaining parameters assuming flat prior probabilities. 
The figures show strong degeneracies as apparent from equa- 
tion ( |42[ ). These degeneracies represent the origin of the 
weakened constraints on ay and a_L, which arise because 
the transmission fluctuations introduce scale and angular 
dependencies into the power-spectrum that mimic those in- 
troduced by an incorrect choice off cosmology. In particu- 
lar, Figure |9] shows that the correlation between an and a± 
nearly disappears. This is because by marginalising over , 
the redshift space distortions are now not perfectly known 
(as was assumed for the traditional galaxy redshift survey 
case), so that the Alcock-Paczynski test cannot be used to 
measure DaH. 

We have also computed the contours of likelihood for 
the parameter sets {Cv,b + Cp) and {bCr,b + Cp), in the 



cases of [Cv , bCr ,{b + Cp 
Cp)] = (0,0,2.5) (Figures [To] and 



(0.5,0.1^.5) and [Cv,bCr,{b + 
respectively). These 



11 



can be considered in addition to the constraints (Ci,,Q||), 
{a±,Cv), {a±,bCr), {bCr,a\\), and represent estimates for 
the available constraints on transmission models. For each 
set the likelihood is marginalised over the remaining param- 
eters assuming fiat prior probabilities. The parameters Cv 



and bCr describing the transmission model show very little 
degeneracy with each other. However constraints on Cv are 
degenerate with the power-spectrum amplitude (b+Cp). The 
constraints on these parameters are similar in magnitude for 
the two fiducial models considered. Without detailed prior 
knowledge of the cosmology (i.e. no prior on a\\ or a±), a 
Lya survey like HETDEX could determine uncertainties in 
Cv, bCr and {b + Cp) of ACv ~ ±0.04, AbCr ~ ±0.2 and 
A(fe ± Cp) ~ ±0.02. We discuss constraints on transmission 
models in more detail in § [6] 



5.7 Cosmological constraints including prior 
probabilities for C„ and bCr- 

We have shown that fluctuations in Lya transmission de- 
crease the precision with which the angular diameter dis- 
tance and Hubble parameter can be measured at z ~ 2.5, 
relative to measurements from a traditional galaxy redshift 
survey. Assuming no prior knowledge of Cv and bCr, the 
decrease is found to be a factor of 1.5-2. On the other hand, 
if the parameters describing the Lya transmission fluctua- 
tions are separately constrained, then the degeneracies seen 



in Figures 10 and 11 11 imply that the parameters ay and a± 
will be measured with increased precision. In this section we 
investigate the degree of prior knowledge regarding trans- 
mission of Lya flux (as parameterised by Cv and bCr) that 
is required in order to achieve measurements of cosmologi- 
cal parameters with a precision that would be available in a 
traditional galaxy redshift survey. 

Figure [12] shows the constraints on power-spectrum dis- 
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Figure 12. Constraints on power-spectrum distortions achievable via the Alcock-Paczynski test including prior constraints on Ct, 
and bCp. Each panel shows contours of likelihood for the parameter set (ox,a||), in the case of a hya galaxy power-spectrum with 
[Cv, bCr, {b + Cp)] = (0, 0, 2.5). Contours of likelihood for the values of ay and ax are shown (at 61% and 14% of the peak likelihood). 
Prior likelihoods of Cq^ {Cv) = exp {—{Cv — {Cv))^ /2a^) and Cqj. {bCr) = exp {—{bCr — {''C'r))^/2tT^p) were included in the constraints. 
The four columns show constraints assuming CTi, = 1, 0.05, 0.02 and 0.01. For each of these columns, four rows are shown with tr^p = 1, 
0.25, 0.1 and 0.05. The 61% contour for the traditional galaxy redshift survey constraints is plotted in each panel for comparison (thick 
dotted lines). 



tortions that are achievable via the Alcock-Paczynski test 
in cases where there are prior constraints on Cv and bCr- 
Each panel shows contours of likelihood for the parameter 
set (ax, ay), in the case of a Lya galaxy power-spectrum 
with [Cv, bCr, (b + Cp)] = (0, 0, 2.5). To compute these con- 
straints we have assumed prior likelihoods of 

CcACv) = exp{-{Cv~{Cv}f/2ai) and 
CcribCr) = exp(-(6Cr-{&Cr))V2<T?r), (45) 

where (Cv) and (bCr) are the means of Cv and bC-, cor- 
responding to the fiducial model, and av and atr are the 
prior uncertainties in Cv and bCr respectively. Figure [l2] 
is presented as a grid, with prior precision on Cv increasing 



from left to right, and prior precision on bCr increasing from 
top to bottom. The four columns show constraints assuming 
(Jv — 1, 0.05, 0.02 and 0.01. For each of these columns, four 
rows are shown with abr ~ 1, 0.25, 0.1 and 0.05. The 61% 
contour for the traditional galaxy redshift survey constraints 
is repeated in all panels for comparison (thick dotted lines). 

As the prior precision on bCr is increased, the precision 
with which a± is measured increases. As the prior precision 
on Cv is increased, the strength of the correlation between 
ay and a± is increased. This signifies that the parameter Cv 
is degenerate with the Alcock-Paczynski effect. As a result, 
precision in ay requires prior knowledge of both bC'r and 
Cv Prior uncertainties with values smaller than atr ^ 0.05 



© 0000 RAS, MNRAS 000, 000-000 



The clustering of Lya emitters 19 



Table 2. Summary of constraints on «|| and aj_ for the different survey's and analysis methods. The label "galaxy" refers to a traditional 
galaxy redshift survey, whereas "Lya gal." refers to an analysis that includes marginalisation over bCr and Cv The labels "BAO" and 
"AP" refer to constraints based on the BAO scale using equation | |40[ l and the power-spectrum shape via equation ( |42| respectively. 
Where relevant the values of the fiducial model Cv and bC'r, and the corresponding prior constraints cry and (T(,p are listed. 
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prior 
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cosmolo 
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BAO 
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BAO 
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galaxy 
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0.70% 


0.85% 


Lya gal. 


AP 












1.25% 


1.35% 


Lya gal. 


AP 


0.5 


0.1 






1.35% 


1.80% 


Lya gal. 


AP 








1.0 


1.0 


1.20% 


1.35% 


Lya gal. 


AP 








0.05 


0.25 


1.05% 


1.10% 


Lya gal. 


AP 








0.02 


0.1 


0.80% 


1.00% 


Lya gal. 


AP 








0.01 


0.05 


0.75% 


0.90% 



and cr„ < 0.01 provide sufficient precision that fluctuations 
in transmission do not dominate the uncertainties in the 
clustering, in which case a Lya galaxy survey could be used 
to measure cosmological parameters with a precision close 
to that available in a traditional galaxy redshift survey. 



5.8 Comparison with previous HETDEX forecasts 



Shoji et al. ( 2009 1 have presented forecasts for a galaxy red- 



shift survey with parameters corresponding to HETDEX. 
Comparison with the results of their study serves both as 
a check of our analysis, and illustrates the relationship be- 
tween the astrophysical parameters introduced through Cp, 
Cr and Cv, and the cosmological parameters a« and a±. 



The results from the analysis of Shoji et al. ( 2009 1 are listed 



in their Table 1, and assume a value of bias 6 = 2.5, making 
the constraints directly comparable to this paper. Firstly, 
our case of a traditional galaxy power-spectrum should be 
compared to the constraints on a|| and a± obtained where 
constraints are marginalised only over the power-spectrum 
amplitude. [Shoji et ah] ( [2009| find Aay = 0.78% and 
Aa± = 0.88% in this case, which corresponds well to our 
values of 0.85% and 0.7%. 

In order to compare our results for a Lya selected galaxy 



redshift survey to the work of Shoji et al. (20091, we first 



consider the case of the limit where fcA <C 1, for which our 
power spectrum model is 

PLy<.(fc, = P{k) [{b + Cp + bCr) + (1 - a)/^'] ' • 

Defining the power-spectrum with a new bias parameter 
such that B = {b + Cp + bCr) and the |Kaiser| ( |1987[ ) factor 
such that F = {1 ~ Cv),f, we obtain 



Pi^yc,{k,fi) = P{k) [B + Ffl 



21 2 



(46) 



which is identical to the usual form. The modifications to the 
bias and Kaiser factor should therefore not affect the Alcock- 
Pacynski effect. However the addition of transmission fiuc- 
tuations means that constraints must be marginalised over 
the Kaiser factor, which changes the shape of contours. 



examples of constraints are shown that include the cases 
of marginalisation over amplitude, and also of marginali- 



sation over both amplitude and the Kaiser (19871 factor. 
Marginalising over the |Kaiser| ( |1987| factor introduces ad- 
ditional uncertainty, with constraints of Aay = 1.13% and 
Aa± — 1.10% in this case. These values should be compared 
to our constraints of Aay — 1.1% and Aa^ = 0.9% which 
are obtained with a broad prior on Cv, but tight constraints 



on C-y (see Figure 12 1. Finally, including the scale depen- 



dent ionizing back-ground term K{k) represents a similar 
effect to that of the primordial spectral index and its run- 
ning. Shoji et al. ( 2009 1 presented constraints that included 
marginalising over the amplitude, [Kaiser] ( |1987[ ) factor, and 
primordial power-spectrum shape, finding Aay = 1.36% and 
Aax — 1.23%. These values are again very similar to our 
constraints without priors on Cp, Cr or Cv of Aay — 1.35% 
and Aax = 1.25% 

Thus, the level at which transmission fluctuations will 
influence measurements of cosmological distance are com- 
parable to those from marginalising over other cosmological 
parameters. However interpretations of measured quantities 
like the power-spectrum shape or the growth function / will 
need to account for the astrophysical effects of Lya trans- 
mission. For example, the possibility of a non-zero Cv will 
complicate interpretation of the extracted value of the red- 
shift distortion factor, because the measurement will be of 
(1 — Cv)f, rather than of /. This degeneracy will make it 
very difficult to test theories of modifled gravity theory using 
the redshift space distortion (e.g. Blake et al. [[201(5 1. 



This can be seen in Figure 3 of Shoji et al. ( 2009J , where 



6 CONSTRAINTS ON LYa TRANSMISSION 
MODELS 

We have shown that fluctuations in Lya transmission will 
result in modiflcation of the scale and angular dependance 
of the power-spectrum measured from a large Lya galaxy 
redshift survey like HETDEX. These modifications in the 
power-spectrum potentially reduce the ability of the mea- 
sured power-spectrum to constrain cosmological parameters. 
On the other hand, our study also shows that the param- 
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Figure 13. Constraints on ttie parameters (b + Cp), Cr and Cy based on power-spectrum distortions. Uncertainties on the cosmology are 
included via the Alcock-Paczynski test. The left, central and right panels show contours of likelihood for the parameter sets {Cv, bCr), 
{Cy, b + Cp) and (bCr, b + Cp), in the case of a Lyo galaxy power-spectrum with [Cy, bCr, {b + Cp)] = (0, 0, 2.5). Contours of likelihood 
for the values of a|| and a± are shown (at 61% and 14% of the peak likelihood). We have included prior likelihoods of Ca{Q:\\) = 
exp (— /2Aa|) and Caj^{ct±) = oxp (— a^/2Aa-^) to represent different precisions of knowledge of the cosmology at z=2.5. In the 
upper and lower panels Aa = Aa = 0.05 and Aa = Aa = 0.01 respectively. 



eters describing the modifications to power-spectrum are 
measured as part of the fitting process. Before concluding, 
we therefore calculate the available precision on the param- 
eters b + Cp, bCr and Cv, whose values provide insight into 
the level of suppression of Lya flux, and the astrophysics of 
the interaction between Lyo emitting galaxies and the IGM 

(§|4}- 



In Figure [13] we present constraints on the parameter 
sets (C„,6Cr), {Cv,b + Cp) and {bCr,b + Cp), in the case 
of a Lya galaxy power-spectrum with [Cv, bCr, {b + Cp)] — 
(0, 0, 2.5). The numerical values of the measured uncertain- 
ties are listed in Table [S] These constraints assume that 
the cosmological model is measured from other sources. The 
uncertainties in the cosmology are included in our analy- 



sis via the Alcock-Paczynski test (equation 421, which for 
this application provides a measure of the uncertainty in 
the power-spectrum shape through the dilation parameters 
Q|| and Q_L. An exception is that the uncertainty in the mass 
power-spectrum amplitude (proportional to the normalisa- 
tion of the primordial power-spectrum, erg) which is degen- 
erate with (6 + Cp). Our analysis assumes that the shape 
of the primordial power-spectrum and the value of / are 
known. 



Using the power-spectrum sensitivity specified in equa- 



tion (39 1, we construct likelihoods 

'PLy(fe,M,P3 - Ply{k,fJ-,Po) 



\nCip) 



E 



APLy(fc,M) 



+ In + hi £c„ + In Ccr 



(47) 



where the sum is over bins of k and fi. We assume flat 
prior probabilities for the transmission dependent param- 
eters (i.e. Ccp ~ -Cc^ ~ jCcr ~ const). To quantify the 
uncertainty we have included prior likelihoods of jC^^^ {a\\) — 
exp (— Qf|/2(T^ ll) and £,a±{a±) = exp {—a\/2a^ j_). In the 
upper and lower panels of Figure[T3]we assume uncertainties 

o'a,± = 0.01 respectively. 



of cr„ 



0.05 and a„ 



Given cosmological uncertainties aa,\\ = Ca.x = 0.05, 
we flnd that the values of b + Cp, bCr and Cv could be con- 
strained with precisions of A(fe + Cp) ~ ±0.02, A(6Cr) ~ 
±0.25 and ACv ~ ±0.04. For next generation cosmological 
constraints with Aan — Aa± = 0.01 we find smaller er- 
rors on Lyof clustering parameters of A(6 + Cp) ~ ±0.015, 
A(6Cr) ~ ±0.15 and AC„ ~ ±0.02. These errors should be 
compared to the predicted values in Table [l] and with the 
results of Figure |4] Since the measurement of non-zero val- 
ues of the parameters b + Cp, bCr or Cv indicates that the 
Lya line is partially absorbed in the IGM (i.e. F > and 
r > in the analytic model), the available constraints in a 
survey like HETDEX could easily measure the presence of 
Lya absorption in the IGM. 
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Table 3. Summary of constraints on {b + Cp) bC'r and Cv The values of the fiducial model (fe + Cp), bC'r and Cv, and the prior 
constraints 0|| and a±_ are listed. 



fiducial model 


prior 


constraint 


Lya transmission constraints 


bCr Cv 






A{b + Cp) A(feCr) A(C„) 





0.05 


0.05 


0.017 0.19 0.040 





0.01 


0.01 


0.013 0.15 0.025 



Of the available constraints, the parameter (6+Cp) can- 
not be used to constrain transmission models despite the 
high precision with which it wiU be determined. This is be- 
cause of the degeneracy between Cp and the unknown galaxy 
bias b. Our models predict small values of bCr (see § 3.5 and 
§[4|, except in very special cases, and so the available preci- 
sion will not provide useful constraints. However the preci- 
sion ACv ~ ±0.04 is small compared with expected values 
of Cv for a range of scenarios (§ [4] see also Zheng et al 
20111. In particular, the HETDEX survey could distinguish 



between infall and outflow dominated models of IGM trans- 
mission, for which the parameters range from 0.05 to 0.7. 
Measurement of this term would directly determine the ex- 
tent to which the IGM impacts the observed flux, providing 
critical information on the intrinsic Lya luminosity, and the 
presence of outflows. 



7 SUMMARY AND CONCLUSION 

Wide-field searches are now finding Lya emitters in large 
numbers, and these galaxies contribute greatly to our under- 
standing of the star-formation history, and of galaxy forma- 
tion. In the near future, very large surveys of Lya emitting 
galaxies will provide precision measurements of large scale 
clustering at 2: ~ 2.5, and allow measurement of cosmolog- 
ical parameters at this previously unexplored epoch. How- 
ever, to realise this goal it is crucial to understand the con- 
tribution to the observed clustering amplitude from fiuctua- 
tions in inter-galactic absorption of intrinsic Lya fiux. In this 
paper we have shown that the environmental dependence of 
Lya absorption can lead to significant non-gravitational fea- 
tures in the redshift space power-spectrum of Lya galaxies, 
under a range of different physical scenarios. We have de- 
rived a physically motivated fitting formula that relates the 
scale and direction dependent Lya power-spectrum to the 
mass power-spectrum [f (fc)] 



X 6 1 + Cr 



arctan 



kX 



+ Cp + (1 - Cv)ffi' 



(48) 



This formula can be used in the power-spectrum analyses of 
a galaxy redshift survey to account for the environmental de- 
pendence of Lya absorption, which includes fluctuations in 
density, ionising background and velocity gradient (parame- 
terised by Cp, Cr and Cv respectively). We have presented a 
simple analytic model to calculate the values of these param- 
eters. Our calculations imply that standard Lya absorption 
scenarios will yield values for Cp, Cr and Cv that are at the 
tens of percent level, indicating that fluctuations in absorp- 
tion of the Lya line in the IGM can lead to modifications 



of the power-spectrum of Lya selected galaxies that are of 
order unity. 

While our analytic model is useful for investigating the 
qualitative dependencies of clustering in Lya selected galax- 
ies, more detailed analyses are required to quantitatively 
predict the values of the constants Cp, Cr and Cv To quan- 
tify the expected effect of fiuctuations in Lya absorption on 
the observed power-spectrum we have therefore employed 
previously published models of Lya radiative transfer. These 
models explore the combined effects of local star-formation 
and IGM infall, and galactic wind driven outflows on the 
transmission of the Lya line through the circum-galactic 
IGM. We find that an infall dominated model for Lya trans- 
mission predicts significant contributions (of order unity) to 
the observed power-spectrum. On the other hand, an out- 
flow dominated model for Lya transmission predicts contri- 
butions to the observed power-spectrum that are an order 
of magnitude smaller, at the level of ~ 5 — 10%. 

We have shown that the expected non-zero values of Cp, 
Cr and Cv will complicate attempts to use the clustering of 
Lya emitters to constrain cosmological parameters in very 
large scale surveys. To quantify the influence of Lya absorp- 
tion on the ability of a large Lya galaxy survey to constrain 
cosmological parameters, we have applied our modifled red- 
shift space power-spectrum to a survey with parameters cor- 
responding to the planned HETDEX. We considered both 
cosmological constraints obtained from the BAO scale, and 
from the the full shape of the power-spectrum as measured 
by the Alcock-Paczynski effect. To base-line our study we 
also consider the case of a traditional galaxy redshift survey 
where the probability of galaxy selection is not a function of 
environment. Our analysis shows that a survey with the pa- 
rameters of HETDEX could measure the BAO scale along, 
and transverse to the line-of-sight with precisions of ~ 3%, 
and ~ 2% respectively. We find that this precision is un- 
affected by modifications to the power-spectrum that arise 
from fiuctuations in Lya absorption. This finding is consis- 
tent with previous studies which have shown that sources of 
scale dependent bias can be removed in order to correctly 
recover the BAO scale. 

For a traditional galaxy redshift survey with the vol- 
ume and galaxy density of HETDEX, much tighter con- 
straints are available through consideration of the full power- 
spectrum shape ( Shoji et al.|2009) . In this case we find that 
find that the Alcock-Paczynski effect could be used to con- 
strain both the line-of-sight and transverse directions at bet- 
ter than the 1% level for a traditional galaxy redshift survey 
in which the shape of the primordial power-spectrum and the 
growth function / are known. However, our analysis shows 
that the dependence of observed Lya flux on velocity gradi- 
ent and ionising background has the potential to compromise 
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the cosmological information available from the full power- 
spectrum shape measured in a Lya selected galaxy redshift 
survey. In a scenario where there is no prior knowledge of 
the details of Lya absorption in the IGM, we find the pre- 
cision of line-of-sight and transverse distance measurements 
the HETDEX would be decreased by a factor of 1.5-2, from 
~0.8% in the case of a traditional galaxy redshift survey, to 
~ 1.3%— 1.7%. The weakened constraints on a and a± arise 
because Lya transmission fluctuations introduce scale and 
angular dependencies into the power-spectrum that are de- 
generate with those of an incorrect cosmology. In particular, 
the effect of a fluctuating ionizing background on the shape 
of the observed power-spectrum of Lya selected galaxies is 
similar to that of an uncertainty in the shape of the primor- 
dial power-spectrum, while fluctuations in velocity gradient 
have an effect that is similar to that of redshift space dis- 
tortions. 

We also investigated the precision with which modelling 
of the Lya radiative transfer must be understood in order 
for HETDEX to achieve a goal of better than 1% distance 
measurements based on Lya galaxy clustering. We find that 
as the prior precision on Cr is increased, the accuracy with 
which a± is measured also increases. As the prior precision 
on Cv is increased, the strength of the correlation between 
ay and a± is increased. However increased precision in ay re- 
quires prior knowledge of both bCr and C„ . Prior uncertain- 
ties with values smaller than ACr < 0.05 and ACv < 0.01 
provide sufficient accuracy that fluctuations in transmission 
do not dominate the uncertainties in the clustering. In such 
cases a Lya galaxy survey could be used to measure cosmo- 
logical parameters with a precision comparable to a tradi- 
tional galaxy redshift survey of equivalent volume and num- 
ber density. On the other hand, these uncertainties are at a 
level below the accuracy with which Cp, Cr and C'v can be 
reliably predicted based on current theoretical understand- 
ing. 

We have turned the above analysis around, and assumed 
instead that the cosmology is known a-priori from other 
sources. The clustering of Lya emitters can then be used 
to measure the impact of the IGM on observed Lya lines, 
and to infer the properties of the Lya transmission model. 
Our models predict small values of Cr , and so the precision 
available from a survey like HETDEX (ACr ~ ±0.15-0.25) 
will not provide useful constraints. However the precision 
ACi, ~ ±0.02 — 0.04 would be small compared with the ex- 
pected values for a range of scenarios (§ |4|. For example, 
clustering of Lya galaxies could distinguish between infall 
and outflow dominated models of IGM transmission. Mea- 
surement of this term in a survey like HETDEX would there- 
fore directly determine the extent to which the IGM impacts 
the observed flux, providing information on the kinematics 
of the cold gas through which the Lya photons are scatter- 
ing. 

The nature of Lya emitters, their role in galaxy forma- 
tion, and their utility as probes of cosmology, and of the 
state of the intergalactic and interstellar media are impor- 
tant topics to which large ongoing observational programs 
promise to make signiflcant contributions. In this paper we 
have shown that power-spectrum measurements from a very 
large survey of Lya selected galaxies could be used to study 
the relationship between the observed Lya flux, and the as- 
trophysics of the galaxy-IGM connection. However we also 



show that Lya transmission fluctuations decrease the cos- 
mological precision of a galaxy redshift survey by a factor of 
1.5-2 relative to the best case precision available in a tradi- 
tional galaxy redshift survey. Realising the full cosmological 
cosmological potential of a survey like HETDEX will there- 
fore require a much more detailed theoretical understanding 
of the astrophysics that determines the relationship between 
the observed and intrinsic Lya flux. Our study underlines 
the need for continued work to understand the role of ra- 
diative transfer through the inter-stellar and inter-galactic 
media in determining the observed properties of Lya emis- 
sion. 

Acknowledgments We thank Eiichiro Komatsu for 
very helpful insights, suggestions, and discussion of our re- 
sults. JSWB thanks the Max Planck Institute fur Astro- 
physik for their hospitality during this work. The research 
was supported by the Australian Research Council (JSBW). 



REFERENCES 

Ahn S., Lee H., Lee H. M., 2003, MNRAS, 340, 863 
Alcock C., Paczynski B., 1979, Nature, 281, 358 
Barkana R., 2004, MNRAS, 347, 59 
— , 2006, MNRAS, 372, 259 

Barnes L. A., Haehnelt M. G., Tescari E., Viel M., 2011, 
ArXiv e-prints 

Blake C., Brough S., CoUess M., Couch W., Groom S., 
Davis T., Drinkwater M. J., Forster K., Glazebrook K., 
Jelliffe B., Jurek R. J., Li I., Madore B., Martin C, Pimb- 
blet K., Poole G. B., Pracy M., Sharp R., Wisnioski E., 
Woods D., Wyder T., 2010, MNRAS, 406, 803 

Blanc G. A., Adams J., Gebhardt K., HiU G. J., Drory 
N., Hao L., Bender R., GiarduUo R., Finkelstein S. L., 
Gawiser E., Gronwall C, Hopp U., Jeong D., Kelzenberg 
R., Komatsu E., MacQueen P., Murphy J. D., Roth M. M., 
Schneider D. P., Tufts J., 2010, ArXiv e-prints 

Bolton J. S., Haehnelt M. G., 2007, MNRAS, 382, 325 

Cassata P., Le Fevre O., Garilli B., Maccagni D., Le Brun 
v., Scodeggio M., Tresse L., Ilbert O., Zamorani G., Cuc- 
ciati O., Contini T., Bielby R., MeUier Y., McCracken 
H. J., Polio A., ZanichelU A., Bardelh S., Cappi A., 
Pozzetti L., Vergani D., Zucca E., 2011, Astron. Astro- 
phys., 525, A143± 

Cole S., Percival W. J., Peacock J. A., Norberg P., Baugh 
C. M., Frenk C. S., Baldry I., Bland-Hawthorn J., Bridges 
T., Cannon R., CoUess M., CoUins C, Couch W., Cross 
N. J. G., Dalton G., Eke V. R., De Propris R., Driver 
S. P., Efstathiou G., Ellis R. S., Glazebrook K., Jackson 
C, Jenkins A., Lahav O., Lewis I., Lumsden S., Maddox 
S., Madgwick D., Peterson B. A., Sutherland W., Taylor 
K., 2005, MNRAS, 362, 505 

Corasaniti P.-S., Huterer D., Melchiorri A., 2007, Phys. 
Rev. D, D75, 062001 

Dayal P., Ferrara A., Saro A., Salvaterra R., Borgani S., 
Tornatore L., 2009, MNRAS, 400, 2000 

Dessauges-Zavadsky M., D'Odorico S., Schaerer D., 
Modigliani A., Tapken C, Vernet J., 2010, Astron. As- 
trophys., 510, A26± 

Dijkstra M., Haiman Z., Spaans M., 2006, ApJ, 649, 14 

Dijkstra M., Lidz A., Wyithe J. S. B., 2007a, MNRAS, 377, 
1175 



© 0000 RAS, MNRAS 000, 000-000 



The clustering of Lya emitters 23 



Dijkstra M., Mcsinger A., Wyithe S., 2011, ArXiv e-prints 
Dijkstra M., Wyithe J. S. B., 2010, MNRAS, 408, 352 
Dijkstra M., Wyithe J. S. B., Haiman Z., 2007b, MNRAS, 
379, 253 

Eisenstein D., 2002, in Astronomical Society of the Pacific 
Conference Series, Vol. 280, Next Generation Wide-Field 
Multi-Object Spectroscopy, Brown M. J. I., Dey A., eds., 
pp. 35-+ 

Eisenstein D. J., Hu W., 1999, ApJ, 511, 5 
Eisenstein D. J., Hu W., Tegmark M., 1998, ApJL, 504, 
L57+ 

Eisenstein D. J., Seo H., White M., 2007, ApJ, 664, 660 
Eisenstein D. J., Zehavi I., Hogg D. W., Scoccimarro R., 
Blanton M. R., Nichol R. C, Scranton R., Seo H.-J., 
Tegmark M., Zheng Z., Anderson S. P., Annis ,J., Bah- 
call N., Brinkmann J., Buries S., Castarider P. J., Con- 
nolly A., Csabai I., Doi M., Pukugita M., Prieman ,J. A., 
Glazebrook K., Gunn J. E., Hendry J. S., Hennessy G., 
Ivezic Z., Kent S., Knapp G. R., Lin H., Loh Y.-S., Lup- 
ton R. H., Margon B., McKay T. A., Meiksin A., Munn 
J. A., Pope A., Richmond M. W., Schlegel D., Schneider 
D. P., Shimasaku K., Stoughton C, Strauss M. A., Sub- 
baRao M., Szalay A. S., Szapudi L, Tucker D. L., Yanny 

B. , York D. G., 2005, ApJ, 633, 560 
Faucher-Giguere C, Lidz A., Hernquist L., Zaldarriaga M., 

2008, ApJ, 688, 85 
Gawiscr P., Prancke H., Lai K., Schawinski K., Gronwall 

C, CiarduUo R., Quadri R., Orsi A., Barricntos L. P., 
Blanc G. A., Fazio G., Peldmeier J. J., Huang J., Infante 
L., Lira P., Padilla N., Taylor E. N., Treister E., Urry 

C. M., van Dokkum P. G., Virani S. N., 2007, ApJ, 671, 
278 

Gaztanaga E., Cabre A., Hui L., 2008, ArXiv e-prints, 
0807.3551 

Giannantonio T., Scranton R., Crittenden R. G., Nichol 
R. C, Boughn S. P., Myers A. D., Richards G. T., 2008, 

ArXiv e-prints, 0801.4380 
Glazebrook K., Blake C, 2005, ApJ, 631, 1 
Guaita L., Gawiser P., Padilla N., Prancke H., Bond N. A., 
Gronwall C., CiarduUo R., Feldmeier J. J., Sinawa S., 
Blanc G. A., Virani S., 2010, ApJ, 714, 255 
Haiman Z., Spaans M., 1999, ApJ, 518, 138 
Hansen M., Oh S. P., 2006, MNRAS, 367, 979 
Heckman T. M., Borthakur S., Overzier R., Kauffmann G., 
Basu-Zych A., Leithercr C, Scmbach K., Martin D. C, 
Rich R. M., Schiminovich D., Scibcrt M., 2011, ApJ, 730, 
5 

Hill G. J., Gebhardt K., Komatsu P., Drory N., MacQueen 
P. J., Adams J., Blanc G. A., Koehler R., Rafal M., Roth 

M. M., Kelz A., Gronwall C, CiarduUo R., Schneider 

D. P., 2008, in Astronomical Society of the Pacific Confer- 
ence Series, Vol. 399, Astronomical Society of the Pacific 
Conference Series, T. Kodama, T. Yamada, & K. Aoki, 
ed., pp. 115 — h 

Hill G. J., Gebhardt K., Komatsu P., MacQueen P. J., 
2004, in American Institute of Physics Conference Series, 
Vol. 743, The New Cosmology: Conference on Strings and 
Cosmology, R. E. Allen, D. V. Nanopoulos, & C. N. Pope, 
ed., pp. 224-233 

Hu E. M., Cowie L. L., Capak P., McMahon R. G., 
Hayashino T., Komiyama Y., 2004, AJ, 127, 563 

Hui L., Gnedin N. Y., 1997, MNRAS, 292, 27 



Ihev I. T., Shapiro P. R., McDonald P., Mellema G., Pen 
U., 2008, MNRAS, 391, 63 

lye M., Ota K., Kashikawa N., Purusawa H., Hashimoto T., 
Hattori T., Matsuda Y., Morokuma T., Ouchi M., Shi- 
masaku K., 2006, Nature, 443, 186 

Kaiser N., 1987, MNRAS, 227, 1 

Kashikawa N., Shimasaku K., Malkan M. A., Doi M., Mat- 
suda Y., Ouchi M., Taniguchi Y., Ly C, Nagao T., lye 
M., Motohara K., Murayama T., Murozono K., Nariai K., 
Ohta K., Okamura S., Sasaki T., Shioya Y., Umemura M., 
2006, ApJ, 648, 7 

Komatsu P., Smith K. M., Dunkley J., Bennett C. L., Gold 

B. , Hinshaw G., Jarosik N., Larson D., Nolta M. R., Page 
L., Spergel D. N., Halpern M., Hill R. S., Kogut A., Limon 
M., Meyer S. S., Odegard N., Tucker G. S., Weiland J. L., 
WoUack E., Wright E. L., 2011, ApJS, 192, 18 

Kovac K., Somerville R. S., Rhoads J. E., Malhotra S., 

Wang J., 2007, ApJ, 668, 15 
Laursen P., Sommer-Larsen J., Andersen A. C., 2009, ApJ, 

704, 1640 

Laursen P., Sommer-Larsen J., Razoumov A. O., 2011, 
ApJ, 728, 52 

Lehnert M. D., Nesvadba N. P. H., Cuby J., Swinbank 

A. M., Morris S., Clement B., Pvans C. J., Bremer M. N., 

Basa S., 2010, Nature, 467, 940 
Lidz A., Paucher-Giguere C., DaU'AgUo A., McQuinn M., 

Fechner C, Zaldarriaga M., Hernquist L., Dutta S., 2010, 

ApJ, 718, 199 
Madau P., 1995, ApJ, 441, 18 
Malhotra S., Rhoads J. P., 2004, ApJL, 617, L5 
Mao Y., Tegmark M., McQuinn M., Zaldarriaga M., Zahn 

O., 2008, Phys. Rev. D, 78, 023529 
Mas-Hesse J. M., Kunth D., Tenorio-Tagle G., Lcithcrer 

C, Terlevich R. J., Terlevich P., 2003, ApJ, 598, 858 
McDonald P., Eisenstein D. J., 2007, Phys. Rev. D, 76, 

063009 

McQuinn M., Hernquist L., Zaldarriaga M., Dutta S., 2007, 

MNRAS, 381, 75 
Mesinger A., Purlanetto S. R., 2008a, MNRAS, 385, 1348 
— , 2008b, MNRAS, 386, 1990 

Miralda-Escude J., Haehnelt M., Rees M. J., 2000, ApJ, 

530, 1 

Mo H. J., White S. D. M., 1996, MNRAS, 282, 347 

Morales M. P., Wyithe J. S. B., 2010, Ann. Rev. Astron. 
Astrophys., 48, 127 

Nagamine K., Ouchi M., Springel V., Hernquist L., 2010, 
PASJ, 62, 1455 

Neufeld D. A., 1991, ApJL, 370, L85 

Okumura T., Matsubara T., Eisenstein D. J., Kayo I., Hik- 
age C, Szalay A. S., Schneider D. P., 2008, ApJ, 676, 889 

Orsi A., Laeey C. G., Baugh C. M., Infante L., 2008, MN- 
RAS, 391, 1589 

Osterbrock D. P., 1989, Astrophysics of gaseous nebulae 
and active galactic nuclei, Osterbrock, D. E., ed. 

Ouchi M., Shimasaku K., Akiyama M., Simpson C, Saito 
T., Ueda Y., Purusawa H., Sekiguchi K., Yamada T., Ko- 
dama T., Kashikawa N., Okamura S., lye M., Takata T., 
Yoshida M., Yoshida M., 2008, ApJS, 176, 301 

Ouchi M., Shimasaku K., FYirusawa H., Saito T., Yoshida 
M., Akiyama M., Ono Y., Yamada T., Ota K., Kashikawa 
N., lye M., Kodama T., Okamura S., Simpson C, Yoshida 
M., 2010, ApJ, 723, 869 



© 0000 RAS, MNRAS 000, 000-000 



24 Wyithe & Dijkstra 



Peacock J. A., Colo S., Norbcrg P., Baugh C. M., Bland- 
Hawthorn J., Bridges T., Cannon R. D., CoUess M., 
Collins C, Couch W., Dalton G., Deeley K., De Propris 
R., Driver S. P., Efstathiou G., Ellis R. S., Prank C. S., 
Glazebrook K., Jackson C, Lahav O., Lewis I., Lumsden 
S., Maddox S., Pcrcival W. J., Peterson B. A., Price I., 
Sutherland W., Taylor K., 2001, Nature, 410, 169 

Percival W. J., Cole S., Eisenstein D. J., Nichol R. C, 
Peacock J. A., Pope A. C, Szalay A. S., 2007, MNRAS, 
381, 1053 

Press W. H., Schechter P., 1974, ApJ, 187, 425 
Reid B. A., Percival W. J., Eisenstein D. J., Verde L., 
Spergel D. N., Skibba R. A., Bahcall N. A., Budavari 
T., Fricman J. A., Fukugita M., Gott J. R., Gunn J. E., 
Ivozic Z., Knapp G. R., Kron R. G., Lupton R. H., McKay 
T. A., Meiksin A., Nichol R. C, Pope A. C, Schlegel 
D. J., Schneider D. P., Stoughton C, Strauss M. A., Sza- 
lay A. S., Tegmark M., Vogeley M. S., Weinberg D. H., 
York D. G., Zchavi I., 2010, MNRAS, 404, 60 
Rhook K. J., Geil P. M., Wyithe J. S. B., 2009, MNRAS, 
392, 1388 

Rybicki G. B., Lightman A. P., 1979, Radiative processes 

in astrophysics, Rybicki, G. B. & Lightman, A. P., ed. 
Santos M. R., 2004, MNRAS, 349, 1137 
Schaerer D., 2003, Astron. Astrophys., 397, 527 
Sec H., Eisenstein D. J., 2007, ApJ, 665, 14 
Sec H., Siegel E. R., Eisenstein D. J., White M., 2008, ApJ, 
686, 13 

Seo H.-J., Eisenstein D. J., 2005, ApJ, 633, 575 

Shapley A. E., Steidel C. C, Pettini M., Adelberger K. L., 

Erb D. K., 2006, ApJ, 651, 688 
Sheth R. K., Mo H. J., Tormcn G., 2001, MNRAS, 323, 1 
Shimasaku K., Kashikawa N., Doi M., Ly C, Malkan M. A., 
Matsuda Y., Ouchi M., Hayashino T., lye M., Motohara 
K., Murayama T., Nagao T., Ohta K., Okamura S., Sasaki 
T., Shioya Y., Taniguchi Y., 2006, PASJ, 58, 313 
Shimizu I., Yoshida N., Okamoto T., 2011, ArXiv e-prints 
Shoji M., Jeong D., Komatsu E., 2009, ApJ, 693, 1404 
Steidel C. C, Bogosavljcvic M., Shapley A. E., KoUmeier 
J. A., Reddy N. A., Erb D. K., Pettini M., 2011, ArXiv 
e-prints 

Steidel C. C, Erb D. K., Shapley A. E., Pettini M., Rcddy 
N., Bogosavljevic M., Rudie G. C, Rakic O., 2010, ApJ, 
717, 289 

Tinker J. L., Weinberg D. H., Zheng Z., 2006, MNRAS, 
368, 85 

Vanzella E., Grazian A., Hayes M., Pentericci L., Schaerer 
D., Dickinson M., Cristiani S., Giavalisco M., Verhamme 
A., Nonino M., Rosati P., 2010, Astron. Astrophys., 513, 
A20+ 

Verhamme A., Schaerer D., Atek H., Tapken C, 2008, As- 
tron. Astrophys., 191, 89 
Verhamme A., Schaerer D., Maselli A., 2006, Astron. As- 
trophys., 460, 397 
Wyithe J. S. B., Loeb A., 2005, ApJ, 625, 1 
Yajima H., Choi J., Nagamine K., 2010, ArXiv e-prints 
Zheng Z., Cen R., Trac H., Miralda-Escude J., 2010, ApJ, 

716, 574 
— , 2011, ApJ, 726, 38 



